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Introduction 



Finite size effects play an important role in modern Statistical Mechanics 
and Quantum Field Theory. In fact, they give a description of the special 
behaviour of observables of a statistical model like specific heat, magnetic 
susceptibility or correlation length due to the influences of the presence of a 
finite geometry. As an example, consider the specific heat c(T) of some statis- 
tical system defined in a finite domain with size L. Although the specific heat 
of the system in an infinite volume shows a divergence Coo(T) ~ (T — T c )~ a 
at the critical temperature T c , this singularity is rounded when the system 
size is finite. The curves c(T) versus T show a maximum at some value of T 
which is shifted away from T c (see Fig. (OD)). One can take the temperature 
of the finite size maximum as an estimate T C (L) of the true critical temper- 
ature T c (oo) = T c and only in a finite interval around this temperature the 
finite size effects are relevant, out of this interval these finite size effects are 
negligible because the correlation length of the system is not longer compa- 
rable with the size. These phenomena are also found to hold true for many 
condensed matter systems. For example, many experimental mesurements 
have been done in physical systems as thin epitaxial layers and multi-layers 
of magnetic substance on non-magnetic substrates (see [if and references 
therein) . 

An important fact is that critical quantities, as the specific heat, have 
a scaling behaviour, i.e. a variation as functions of L, that is fixed by the 
critical exponents of the system in the infinite size geometry. Moreover, 
the behaviour of the scaling functions (see later) is deeply related to the 
conformal field theory describing the critical points of statistical systems 
(see |2J and [ll, [3] for a review). 

From the point of view of quantum field theories defined in a finite geom- 
etry, interesting phenomena appear as well. Take, for instance, a 1 + 1 dimen- 
sional space-time with periodic boundary conditions in the space direction 
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Figure 1: Specific heat as a function of the temperature for a system in a 
finite volume L. 



and an infinitely extended time direction, in fact, this geometry corresponds 
to a cylinder: there will be some Casimir effect changing, for example, the 
energy of a two body interaction, because two particles can interact in two 
possible directions along the space direction due to the periodic boundary 
conditions (Fig. (0)). 

New radiative corrections to the self-energy of a propagating particle 
should appear as well. While in an infinite space-time only local virtual 
emissions are allowed, on a cylindrical geometry one can also conceive vir- 
tual emissions traveling around the whole cylinder before coming back to the 
bare particle. Although these virtual emissions are exponentially depressed 
with respect to the usual local emissions, if the cylinder circumference L (the 
size of space direction in Fig. (p£)) is small enough, there is a mesurable 
probability that such virtual objects travel around the world and come back 
from the other side, due to the Heisenberg principle. The corrections to the 
self-energy can be estimated. 

Useful applications of such results show up in lattice calculations, for ex- 
ample many progresses have been made in some two dimensional quantum 
field theories sharing phenomena with QCD and other interesting four di- 
mensional theories. In a particular class of two-dimensional quantum field 



Introduction 




Figure 2: Particles a and b interacting in two different directions in a 
1+1 dimensional space-time with finite space size L and periodic 
boundary conditions. 



theories, the so called Integrable quantum field theory (IQFT), there are 
methods to exactly calculate the dependence of physical quantities on the 
finite size L. 

In principle, the behaviour of such systems in finite geometries can be 
known. This will be the main topic of the present thesis. Among all possible 
methods to compute finite size effects, we will propose the Nonlinear Integral 
Equation approach to study the sine- Gordon model in a finite volume with 
Dirichlet boundary conditions. This approach will be developed through a 
lattice regularization of the model in order to describe the ground state, i.e. 
the Casimir energy, and some excited states scaling functions to give predic- 
tions on the spectrum in relation to the space size L. It will be shown that in 
the ultraviolet limit the Nonlinear Integral Equation allows to reproduce the 
correct conformal field theory living in the ultraviolet renormalization group 
flow fixed point of the model and to classify conformal states in relation to 
physical states of sine-Gordon. 

The thesis is organized following the scheme: in Chapter 1 some known 
fundamental facts about critical phenomena and conformal field theories are 
introduced. Particular relevance is reserved to the relations between statisti- 
cal mechanical systems at criticality and conformal field theory, to conformal 
field theory defined in a cylindrical geometry and to conformal theories with 
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Figure 3: Virtual emissions from the particle a. Due to periodic boundary 
conditions in the space direction, they can travel along the whole 
space coming back to the particle from the left. 



boundaries. Chapter 2 contains a general survey of IQFTs and boundary 
IQFTs. The S matrix approach and the finite size scaling in relation to con- 
formal invariant theories are introduced. Chapter 3 is devoted to introduce 
basic material on the sine-Gordon model with Dirichlet boundary conditions 
and the inhomogeneous XXZ spin chain with boundaries. A description of 
double row transfer matrix, its diagonalization and the derivation of Bethe 
Ansatz equations is also provided. Finally, in Chapter 4 the Nonlinear Inte- 
gral Equation for our model is derived while the analysis of the continuum 
theory is proposed in Chapter 5, where the finite size effects of boundary 
sine-Gordon model are also discussed taking into account excited states. 
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Introduction 



Chapter 1 



Scaling and conformal field 
theory 



This Chapter presents a summary of critical phenomena in statistical me- 
chanical systems, conformal field theory and their mathematical and physical 
relations. 



1.1 Critical phenomena 

A statistical mechanical system is said to be critical when the correlation 
length £, i.e. the distance over which the order parameters are statistically 
correlated, diverges: £ — ► oo. This is peculiar of continuous (second order) 
phase transitions, which consist in a change of macroscopic properties of 
the system when some parameters, for example the temperature, are varied. 
The special points in the phases diagram characterizing a second order phase 
transition are called second order critical points. In particular, a diverging 
correlation length makes no more relevant the length scales, i.e. every phe- 
nomenon at length scale p < £ is on average independent of the value of 
p. Then the scale invariance emerges. As a consequence, the microscopic 
details, such as the precise structure of the interactions, are not longer dis- 
cernible. Therefore, the behaviour of the observables in a region close to 
second order critical points, which typically can be described by critical ex- 
ponents, becomes independent of microscopic properties. The independence 
of quantities like critical exponents from microscopic aspects of a system is 
referred as universality. Phenomenologically, the various critical behaviours 
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of systems near their critical points can be organized in universality classes, 
which depend only on the space dimensionality and on underlying simme- 
tries. One would like to be able to define a systematic classification of the 
possible universality classes, a problem in general not yet achieved. For 
two-dimensional systems, however, conformal invariance allows a satisfying 
description of universality classes. Indeed, the use of conformal invariance to 
describe statistical mechanical systems at criticality is motivated by an argu- 
ment, due to Polyakov [|4j, which states that local scale invariant field theories 
are conformally invariant. Hence, universality classes of critical behaviours 
can be identified with conformal field theory (CFT), i.e. a conformally in- 
variant quantum field theory. 



1.1.1 Scale invariance and scaling 

First of all let us introduce the concept of critical exponents. Given a statisti- 
cal mechanical system, one can find several macroscopic "quantities" deeply 
describing the system, such as the temperature or the magnetic field, if any. 
In particular one defines the "reduced quantity", which describes the devia- 
tion of the system from critical value, e.g. for the temperature one defines 
the reduced temperature as 

T — T 
t=~^ (LI) 

J- c 

where T c is the critical temperature, i.e. the temperature of the critical point. 
Given the basic quantity of a statistical system as the partition function 

z = J2 ex p (-^ H ) ( L2 ) 

states ^ ' 

with H the Hamiltonian, averages (X) for several thermodynamic observ- 
ables can be obtained from it as functions X(t) of the reduced temperature 
t. Away from criticality it is known that the functions X(t) follow expo- 
nential law behaviours X(t) ~ \t\ v± , where the critical exponents y± , s are 
defined as 

y± = hm — — — — (1.3) 

y t-»±o In t v ; 
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In most cases y + = y_, the exponents are the same for t > and t < 0. The 
most common critical exponents of statistical mechanical systems are those 
associated to the specific heat C, to the spontaneous magnetization M and 
to the susceptibility x- They are defined as follows 

C~\t\~ a , M~\-tf , x~l*H ( L4 ) 

We are mostly interested in the critical exponent v of the correlation length 
and the large distance behaviour of the two-point correlation function Giv) 
defined usually as 

i ~ \t\' v , G(r) ~ r 2 - d ~ r > at t = (1.5) 

Critical exponents describe phenomenogically the behaviour of statistical sys- 
tems close to criticality. It is not simple in general to give an exaustive de- 
scription of critical exponents, while in two-dimensional critical phenomena 
the conformal invariance of correlation functions at criticality allows one to 
find the exponents exactly. 

In fact the critical exponents can be related to each other by use of the 
scaling hypothesis, which states that the free energy density near the critical 
point is a homogeneous function of some parameters, such as the external 
magnetic field h and the reduced temperature t, called in this framework 
scaling fields. That is there exist parameters at and ah such that 

f(t,h) = \~ d f(\ a %\ a »h) (1.6) 

where X~ d is a dilatation factor describing a scale transformation of param- 
eters t and h, a t and ah are called renormalization group (RG) eigenvalues. 
From this important properties the critical behaviour of scaling fields can 
be derived. To be generic, one can take a scaling field and the associated 
RG eingenvalue u. Consider the scaling relation $ = \ u (p for the free energy 
/, where the critical point corresponds to = $ = 0. If u > 0, repeated 
dilatations will move the system away from criticality, given an arbitrarily 
small initial perturbation <p ^ near the criticality. Such a scaling field is 
called relevant. If u < 0, repeated dilatations will leave the system around 
the critical point. In this case the scaling field is called irrelevant. Finally, if 
u = 0, the scaling field is said marginal. It follows that the behaviour of the 
system is completely determined by the relevant scaling fields. In particular, 
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for the free energy the scaling hypothesis implies the homogeneity relation 



f(t,h) 



d/at f / . \ / y. \ -a-h/at 



6± 1 1 sJ U i i (L7) 



where ^ is the scaling function and the index ± refers t > or t < respec- 
tively. The scaling function is universal, i.e. independent of the irrelevant 
scaling fields. 

Evidently, because the most relevant thermodynamic quantities of a sta- 
tistical system can be obtained from the free energy density as derivatives 
with respect to t and h, the scaling hypothesis gives a tool to calculate all 
critical exponents as functions of a t and a^: 

o d R d-a h 2a h -d 

a = 2 , /3= , 7 = (1.8) 

a t a t a t 

For the scaling of a two-point correlation function a scale transformation 
imposes 

G(r, t, h) = X~ 2x G (j, X at t, X ah h) (1.9) 

where x is called the scaling dimension, then, eliminating the dilatation factor 
fixing X at t = K and taking for simplicity h = 0, one has the homogeneus 
relation 

G(r,t,0) = t 2 -/-G(^ 7 ^ 7 ,^,o) (1.10) 

Close to criticality, G(r) decays exponentially with |r|, so one can identify a 
correlation length £, which shows the scaling 

£~£ 1/at therefore u = l/a t (1.11) 

Finally, considering the limit t — > 0, one gets 

r] = 2x t + 2-d (1.12) 

Observe that, given the two last relations, the critical exponents a , /3 and 7 
can be found by substitution in the formulae ( 11.81 ). Let us mention that in 
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the particular case of t = h = 0, a spatial dilatation (scale transformation) 
r — ► r' = r/A gives the following relation for two-point correlation functions 

G(r) = A- 2a 'G(r/A) (1.13) 

At least heuristically, for a field theoretic interpretation of this, it appears 
natural to assume that under scale transformations any field obeys the scaling 
behaviour 

&(r) = A-*Vi(r/A) (1.14) 

with Xj scaling dimension of the field 0j. The scaling hypothesis and all the 
procedures illustrated can be motivated and rigorously derived by the renor- 
malization group theory. An exhaustive presentation of the renormalization 
group is out of the scope of this work we just present here a brief survey. 



1.2 Renormalization Group 

A quantum field theory is in general not invariant under scale transforma- 
tions, so it is a very important task to determine the behaviour of a theory 
when scale transformations are performed: this corresponds to study the 
renormalization group (RG) of the theory. One-way to investigate the RG 
is to introduce the Space of Actions of Wilson |5|, |6|. It corresponds to sup- 
pose that the action of a given theory A = J d d x£ depends on a set of 
fields and their derivatives and on the coupling constants g = (gi, . . .g n ). 
The central idea is that one can see the variation of the Lagrangean den- 
sity under scale transformations x — > x' = x + xdt as a transformation 
of the couplings g — ► g'. Scale transformations induce a motion on the n- 
dimensional space of coupling constants. Such a space is usually called Space 
of Actions. A trajectory in this space is thus a function g(t) of the scale pa- 
rameter t. Under a scale transformation the Lagrangean density changes as 
£(t + dt) = C{t) + dfj, J M where J M = x v T^ u is the conserved current and T^ v 
is the energy-momentum tensor. If the system is scale invariant the quantity 
J^ is conserved and then T.f = O = 0, of course this situation does not 
appear for a general g(t). Let us introduce the /5-function 

A(g) = ^ (1-15) 

11 
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which describes the variation of coupling constants under a scaling of the 
parameter t. Given the /3-function for a certain initial value of couplings go, 
as t increases there are different behaviours 

• if (3 (go) > 0, then g(t) is increasing in a neighbourhood of go; 

• if (3 (go) < 0, then g(t) is decreasing in a neighbourhood of g . 

There exist special points g* called fixed points where (3 (g*) = and in- 
dipendently of the exact value of go in a given region, the asymptotic value 
of g(t) for t — > ±oo depends on the values of g closest to g*. In particular, 
one can distinguish two types of fixed points: infrared (IR) fixed points, if 
they are reached from go for t — > +oo; ultraviolet (UV) fixed points, if they 
are reached for t — » — oo. Usually the trajectories g(t) are referred as Renor- 
malization Group flows. The control of /3-functions allows one to reconstruct 
the behaviour of the theory around the starting point go and hence a de- 
scription of the tangent space to the space of actions. In the tangent space 
the variation of the Lagrangean can be described by combinations of fields 

*i = ^ (1-16) 

In particular one has that the trace of the energy-momentum tensor is a field 
living on the tangent space 

2*(x) = 0(x) = Zf^ = E»(B) ( L17 ) 

i i 

Therefore 0(x) = iff /9j(g) = V i G N. This means that the fixed points 
of a quantum field theory are scale invariant. This is a very important result. 
In the renormalization theory it is possible to give a description of the 
variation of N point correlation functions along the RG flow as well. What 
one gets is the so-called Callan-Symanzyk equation. We do not illustrate 
the general derivation (for a nice procedure see |7J) and focus on the two- 
dimensional case. Given a set of N fields *4;(x) and their N point corralation 
function (X) = («4i(x). . . ..Ajv(x)), any field changes under a scale transfor- 
mation as 

^■(x) = dt (rfA. + D,) A(x) (1.18) 



dxt 
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where Dj is the classical dimension. One can demonstrate that for a two- 
dimensional quantum field theory the correlation function obeys the following 
equation along the RG flow 



£(^ + r ^>)-I>< 



(X)=0 (1.19) 



where Tj(g)Aj(x) = [Dj + X/iA(s)^) *4j( x ) is the anomalous dimension 
of the field Aj. Moreover, the dimension of the field $i(x) is 



E( 2 ^ 






r(g)$,(x) = > ][28l - -^ $,(x) (1.20) 



and, in particular, T(g)G(x) = 2Q(x). Therefore, considering also the 
Lorentz invariance, it follows that all the components of the energy- momentum 
tensor have the same anomalous dimension: 2. 



1.2.1 Phenomenological renormalization 

Consider now a statistical mechanical system, say a spin chain, defined in 
a finite (one dimensional) lattice with dimension L in which, for simplicity, 
h = 0. According to finite size scaling, for any change of scale, the correlation 
length obeys the homogeneous relation 

C\t,o) = L~ l g (L^%0) (1.21) 

Given two lattice sizes L and L' and temperatures t and t' such that it is 
satisfied the condition 

L - L> [ ] 

one can describe the relation between t and t' as 



t' = t(-) (1.23) 
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These relations can be reinterpreted [8] as a RG mapping from the two tem- 
peratures 

t^t' = K L , b {t) (1.24) 

where b = L/L' is the scale factor. This procedure can be considered exact 
only in the limit L, L' —* oo, but the definition of the RG mapping is con- 
sidered to make sense also for L, V finite and large enough. If it is found a 
fixed point t* (L, L') of the RG flow, i.e. 

mL± = mn) (125) 

one expects that t* (L, L') — > as L, L' — >• oo. The solution can be derived 
in several ways, for example graphically. The important thing is that the 
RG approach (as it has been done in the previous discussion) allows one to 
relate critical (scale invariant) points to fixed points and then the change of 
physical quantities as scales are continuously varied. From another point of 
view, one can use the fact that the function Q{t,h) is universal and obtain 
that 

SM = (ir) (126) 

Therefore if the value £(0,0) is known, one can estimate t* and give an exact 
description of critical exponents. 



1.3 Conformal invar iance 

It has been said that a CFT is a quantum field theory endowed with covari- 
ance properties under conformal transformations (see |9|, [10j), i.e. special 



coordinate transformations in d dimensions which keep the metric invariant 
up to a scale factor 

V ( r ') = Q ( r )9^ ( r ) 
The set of conformal transformations forms a group: the conformal group. In 
d > 2 this group is finite dimensional and it is generated by global rotations 
r — ► r' = Ar, global translations r — ► r' = r+a, global dilatations r — ► r' = br 
and the so called special conformal transformations 
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. r + ar 2 

r 



1 + 2ar + a 2 r 2 

We restrict to the two dimensional case and in particular to the Euclidean 
plane where the conformal transformations are realized by any analytic change 
of coordinates z — > £(z), z — >• £(z). In fact, one can demonstrate |9| that the 
contributions of the complex variables z and z decouple and that they can 
be regarded as independent variables. One can identify a conformal theory 
as a quantum field theory described by the correlation functions of a set of 
local scaling operators A(z, z) with following properties 

• If A(z, z) is local, all derivatives of A are local. The set of operators 
{A(z,z)} is in general infinite; 

• there exists a subset {§(z,z)} C {A(z,z)} of local operators, called 
quasi-primary operators, which transform covariantly under projective 
conformal transformations z — > w(z) = (az + b)/(cz + d) as 

*<-->- en© ~ %m ■ A ' AeR; 

• any local operator can be written as a linear combination of quasi- 
primary operators and their derivatives; 

• the vacuum of the theory is invariant under projective conformal trans- 
formations. 

Under any conformal change of variables, the fields of a CFT obey some 
general relations. In particular primary fields, which are a subset of quasi- 
primary fields, transform according to 

where the real numbers (A, A) are called conformal dimensions of the field 0. 
This relation implies that it is very simple, in principle, to calculate exactly 
the correlation functions among primary fields using the fact that 
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(01 (wx, w x ) ... <p n (w n , w n )) = 

= n(^) (^J (M*l,2l) ...4>n(Zn,Z n )) (1.28) 

For instance, the two point correlation function explicitly reads 

{<Px (z u z x ) 2 (z 2 , z 2 )) = 6 1 2Zx 2 £z 2 xi (1.29) 

where Zy = Zi — Zj, A x = A 2 = A, A x = A 2 = A. The correlation vanishes 
if the conformal dimensions of the two fields are different. It is evident that 
the expression of the two point function allows one to identify the critical 
exponent 77 given in Eq. ( 11.121 ) as 

r] = 2A + 2-d = 2A (1.30) 

In this scheme, i.e. studying the conformal properties of a massless field 
theory living in a scale invariant second order critical point, one finds exactly 
the critical exponents of the system. Similar analysis can be done for three 
and four point functions. 

At this point it is useful to introduce the energy- momentum tensor T^ u . It 
is defined classically considering the variation of the action under coordinates 
transformations, in quantum theory one makes use of Ward identities. It has 
some very important properties: the symmetric property T12 = T21 as a result 
of rotational invariance, and the tracelessness T11+T22 = as a result of scale 
invariance. In addition the energy-momentum tensor is always conserved, i.e. 
dpT^ = 0. In general one introduces complex components 

/ T zz = l(Txx - T 22 - 2iT 12 ) 

\ T- zz = \{T ll -T 22 + 2iTx 2 ) [i - 6i) 

while T zz = T zz = | (T u + T22) = because of tracelessness. The zz compo- 
nent of the energy-momentum tensor is analytic as a consequence of trace- 
lessness and conservation and denoted by T(z). The zz component T{z) is 
antianalytic. T(z) is the generator of infinitesimal conformal transformations 
z — > w = z + e(z) and likewise T(z) for z — > w = z + e(z), in the sense that 
under this change of variables a correlation function changes as 

16 
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(5(0i (Zl,£l). . . <j)n{Zn,Z n )) = 

= y~ f dxe(w)T(w) (0i Oi, Zi) . . . 4> n (z n , z n )) + c.c. 

with an integration contour C encircling all points Z\. . . z n . The change of 
T(z) itself has the form 



dw\ 



-2 



dz J 



c 



T(w)=[— T(z)--{w,z} (1.32) 
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where {w, z} denotes the schwarzian derivative 

a 3 w o /a!n x : 
I w z \ = d*L _ £ J d*_ 

L ' J dw o I dw 

dz \ 9^ 

with the constant term c called central charge. In simple words the energy- 
momentum tensor T(z) transforms under analytical changes of coordinates 
as a primary field with conformal dimensions (2, 0) up to the schwarzian 
anomaly. Assuming that correlation functions are well defined in the complex 
plane with a finite number of possible singularities, there exist short distance 
expansions of products of fields usually called operator product expansion 
(OPE) 

T(w)(j)(z, z) = ^-r H h regular terms (1.33) 

(w — z) w — z 

m / ^ m / x c/2 2T(z) dT(z) , . ,, 

T(w)T(z) = - — '—— + Mrr + ^ + regular terms (1.34) 

(iw — -s;) 4 (w — zy w — z 

that describe the singular behaviour of correlation functions (T(w)<j)(z, z). . .) 
and (T(w)T(z). . .) as z — > w. 

1.3.1 Virasoro Algebra 

The picture given up to now is the one dealing with correlation functions, 
but, as usual in quantum field theory, it is possible also to deal with an 
operator formalism that describes the system by states, i.e. as vectors in a 
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Figure 1.1: Mapping from an infinite space-time strip to the complex 
plane. 



Hilbert space. In CFT it is useful to work with radial quantization in the 
complex plane, that is surfaces of equal time are circles centered at the origin 
and the Hamiltonian is the dilatation operator. It is denoted by the mapping 
z = e Mt-i*)/L^ % _ e 2ir(t+ix)/L w ^ f euclidean time (see Fig. (JTTTT)). 

On any circle there is a description of the system in terms of a Hilbert 
space of states on which field operators act. 

At this point one can expand the energy-momentum tensor (remember 
that it is a conserved current) on its Laurent modes 

+00 
T(z) = J2 z ~ n ~ 2L n (1-35) 

n=— 00 

and, given a field A(z,z), putting this expansion into the OPE provides a 
definition of what L n A(z) is. For instance L A(z) = AA(z), L-iA(z) = 
dA(z), etc. In general one does not expect fields with negative dimensions, 
then, for every field, L n A(z) must vanish for n positive large enough. For 
example, the primary fields have the highest singularity behaving as 1/z 2 , 
they satisfy L n A(z) =0, n > 0. Intuitively the L n are operators acting 
on the space of fields of the theory, then it is tempting to ask oneself what 
the algebra of these operators is. Commuatators among operators L n can be 
calculated considering the products L n L m A and L m L n A. Making use of the 
OPE among energy-momentum tensor and fields and after integrations over 
closed contours via the Cauchy theorem, one has 

[L n , L m ] = {n- m)L n+m + — n (n 2 - l) S n+mfi (1.36) 

This is the celebrated Virasoro Algebra that is an infinite dimensional Lie 
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algebra in which the central charge plays the role of the coefficient of the 
central extension. Observe that there exist similar commutation rules for 
the antiholomorphic part of the theory and, in particular, any commutator 
between operators L n and L m vanishes because the holomorphic and anti- 
holomorphic parts of the fields are decoupled. 

1.3.2 Representation theory of the Virasoro algebra 

A very important use of the Virasoro algebra is to provide one with a natural 
structure to organize all the states, and then all the fields, of a given CFT. 
The representations which one shall consider are the highest weight repre- 
sentations or Verma modules. They are characterized by a real number A. 
A Verma module V C (A) is generated from a highest weight vector of a CFT 
with central charge c denoted |A) satisfying 

L |A) = A|A) , L n |A) = 0, VneN (1.37) 

by the action of the operators L_ n , n > 0. Among these eigenstates there is 
the vacuum |0) and 



V C (A) = Span {L_ fcl , L_ fc2 , . . . L_ fcr } , I < k x <k 2 < . . . < k r (1.38) 

The whole set of Hilbert states is thus organized into products of represen- 
tations of the holomorphic and antiholomorphic Virasoro algebras, for which 
primary fields are highest weights. In compact form this is usually expressed 

as 

n = A/" AiA V c (A) ® V C (A) (1.39) 

A,A 

Such a representation has the important property of being graded for the 
action of the Virasoro generator L . This means that the spectrum of L in 
V C (A) is {A, A + 1, A + 2, . . . }. The numbers A/a, a count the multiplicity 
of each representation in 7i, this implies they must always be non negative 
integers and A/a, A = if a certain representation V C (A) ® V C (A) does not ap- 
pear, but they are not completely fixed by conformal invariance. Constraints 
on them come from other physical requirements such as locality [ill , 12] or 
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modular invariace [2] • The situation is then quite similar to the case of an- 
gular momentum in ordinary quantum mechanics where the space of states 
can be organized in terms of representations of the angular momentum alge- 
bra, of course here we have an algebra with an infinite number of generators. 
Note that the operators L , L±, L_\ form a subalgebra of the Virasoro al- 
gebra and, in particular, L + L is the generator of dilatations: this gives 
a representation of the Hamiltonian operator in term of Virasoro generators 
as H ~ L + L . 

The correspondence between the primary fields and heighest weight vec- 
tors is given by 

|A,) = Ai (O)|O) (1.40) 

All other fields of the theory associated with states L_ n | A) can be obtained 
from primary fields as descendants, in the sense that <p~ n {z) = (L_ n 0) (z) 
with conformal dimensions A + n. They form, together with primary fields, 
the so-called conformal families [0aJ- Any conformal family forms a repre- 
sentation of the Virasoro algebra, because under conformal transformations 
each member of the family is mapped into a representative of the same family. 
Given a module V C (A) it could be interesting to ask oneself whether this 
module is unitary, i.e. characterized by the absence of negative-norm states 
and whether it is reducible or not, in the sense that it contains a singular 
vector (usually called null vector), i.e. a vector satisfying the axioms ( 11.371 ) 
of a highest weight vector. If so, from this vector a highest weight module 
descends, which is a submodule of V c ( A) too. The Verma module can contain 
several such submodules with non trivial intersections. One constructs the 
irreducible representation by quotienting out the submodule(s) of V C (A). Of 
course, one has to define systematically the ranges of c and A where unitary 
irreducible representations of the Virasoro algebra are possible. It has been 
shown in [13] that unitarity is only possible in the following two cases: 

• for c > 1: all representations are unitary and unitarity only implies 
A > 0. However, the number of primary states is infinite. 

• for < c < 1: the following set of V C (A) is unitary 

c= l ? A _ (r(x+l)-sx) 2 -l 

x(x + l) 4x(x + l) 
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with x G Q and r,s G N s.t. 1 < r < s < x. CFTs with central charge and 
conformal dimensions in the set ( jl.411 ) have a finite number of primary states 
and consequentely a finite number of conformal families. They are called 
minimal models. The structure of the Verma modules defined by primary 
fields 4>< r>s ) is encoded in the corresponding Virasoro characters in terms of 
which the partition function is exactly known 



X 



M = T c/24 Tr q Lo = g~ c / 24+A - £) d r , s (n)q n (1.42) 



n=0 



where d r)S (n) is the number of linearly independent states of the representa- 
tion {(j)(r,s)} at level n. The meaning of the variable q will be clear in the 
next sections. 

1.3.3 Correlation functions 

All fields A(z, z) of a CFT obey a closed algebra under OPE 

Mz, z)Aj{0, 0) = J2 B ijMQ, 0) (1.43) 

k 

where the indices i,j, k run over all fields of the theory. Conformal invari- 
ance reduces the problem to the computation of the OPE algebra among 
primary fields, indeed correlators among descendant fields can be reduced to 
correlators among primaries because of Eq. ( J1.33J ). The OPE algebra among 
primary fields reads 

0,(^,^)0,(0,0) = ^ C ^A fc -A ! -A,-A fe -A,-A J [0fc(O)Q)] (1 44) 

k 

where i,j, k count only primaries and [0^(0, 0)] denotes, as usual, the confor- 
mal families with ancestor the primary 0&. If the constants C^ are known, it 
is possible to reduce all the correlators to the two and three point functions. 
The two point function has been written in Eq. ( 11.291 ) and for the three point 
function one has, making use of the projective invariance only, 

(0i (*i, zi) 02 (z 2 , z 2 ) 03 (z 3 , z 3 )) = C^zlfzlfzlfzlfzlfzlf (1.45) 

where 7^ = A^ — Aj — Aj with i ^ j ^ k and i,j,k = 1,2,3. Therefore, 
almost in principle, all correlation functions are known. If the Cf-'s are 
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unknown, constraints come from the requirement of associativity of the OPE 
among fields, which corresponds to impose duality properties to the four 
point functions. 



1.4 CFT on a cylinder 



It has been observed that in two dimensions the conformal transformations 
are represented by analytical transformations, so they are in general infinite 
in number. Among all possible anlytical transformations there is one which 
will be very important in next Chapters. It is the strip geometry defined by 

z ^ w= l n ( z ) (1.46) 

which maps all the plane onto a strip of width L. The negative part of 
the real axis corresponds to both edges of the strip, which has therefore the 
topology of a cylinder (periodic boundary conditions in the space direction). 
Under this particular transformation Eq. ( 11.321 ) can be rewritten as 

T *M = (t)\ t{z) * 1 ~ £) (L47) 

As a consequence, the expectation value of T in the cylinder is always differ- 
ent from zero; for instance, if (T(z)) = 0, 

(T cy iW) cyl = -^ (1.48) 

The physical meaning of the apparence of the central charge c, also known as 
conformal anomaly, is the way a specific system reacts to macroscopic length 
scales introduced, e.g. boundary conditions. It gives a direct mesure of the 
variation of the energy, or free energy, due to finite size effects. Due to the 
correspondence between the free energy of two-dimensional systems and the 
ground state energy of one-dimensional quantum systems, one can realize the 
analogy with Casimir effect in quantum electrodynamics (there H plays the 
role of c) . In the scheme of a CFT on a cylinder, the Hamiltonian reads 
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Another quantity of interest is the two point function of a primary field 
<(>(z,z) with conformal dimensions (A, A) . In order to write it down on the 
cylinder one needs to use the covariance relation ( 11.271 ) of primary fields with 
the mapping ( 11.461 ) and gets 



(<f>{w 1 ,W 1 )<f>(W2,W 2 )) 



2vr 

T 



LA 



4 S mh 7r ^ 1 - U,2j smh 7r ^ 1 " ti?2j 



-2A 



1.50) 



where it has been supposed, for simplicity, A = A. Observe that for \w\ — 1U2I 
and \w\ — W2I much smaller than L the effect of the finite size disappears and 
one recovers the infinite plane result. In the opposite case one can easily see 
that the correlator decays exponentially with a correlation length 

£ = sk (L51) 

The apparence of this correlation length is due to the presence of the size 
scale L. In a more general way, one takes a set of scaling operators (pi with 
scaling dimensions x% . It has been observed that they scale as in ( 11.141 ) and for 
any two point correlation function a correlation length £, can be introduced. 
Therefore, at criticality, making use of conformal invariance, it is simple to 
identify exactly the scaling dimensions of fields (in the special case of scalar 
primary fields: Xi = 2A$). Moreover, with respect to the formulae ( 11.131 ) and 
one can show that for an infinitely long strip of finite width L and 



with periodic boundary conditions the scaling functions of fields 0j at the 
RG fixed points (critical points) are exactly known 

0(0, 0) = Aizxi = 4?rA i . (1.52) 

CFT gives an a posteriori confirmation of the scaling hypothesis and, from 
a physical point of view, what which one is more interested in, i.e. the 
quantitative relations among universal scaling amplitudes, critical exponents 
and correlation functions. 

1.4.1 Modular invariance 

Up to now it has been assumed that the CFTs live in the whole complex plane 
or in the cylinder described in the previous section. It could be useful to study 
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these theories on a torus geometry. The motivation to do that is the pos- 
sibility of extracting constraints on the content of the theories coming from 
the interactions of the holomorphic and antiholomorphic sectors revealed by 
the so called modular transformations. In this scheme the Hamiltonian and 
the momentum operators propagate states along the two different directions 
of the torus and the spectrum is embodied in the partition function. 

A torus may be defined by two linearly independent vectors on the plane 
and identifying points that differ by an integer combination of these vectors. 
On the complex plane these lattice vectors can be represented by two com- 
plex numbers U\ and u^. The relevant parameter is the ratio r = uj\/uj2- 
Regarding the torus as a cylinder of finite length whose ends are glued to- 
gether, one has an expression for Hamiltonian as the one given in Eq. ( 11.491 ). 
i.e. H = (2tt/L)(Lq + Lq — c/12) and for the momentum operator one has 
P = (27ri/L)(Lo — Lq). Choosing then u\ real and equal to L, it is possible 
to write the partition function as 

Z(t) = Tr ^o-c/24^ -c/24\ ( L53 ) 

where 

q = e 2niT , q = e~ 2nif (1.54) 

For a CFT it is necessary to have the partition function Z{r) invariant under 
the modular group SL (2,Z) generated by 

T : r^T + 1 , S : T->— (1.55) 

T 

This invariance imposes constraints on the operator content of the theory. For 
example, given the Virasoro algebra character functions Xa{q), the partition 



function is [2, 14 

z(q,q) = ^Ca^axaWxa® ( L56 ) 

A, A 

with A/"a A non negative integers characterizing the operator content of the 
theory, therefore its class of universality. It is possible to demonstrate that 
the requirement of modular invariance fixes the numbers A/a, a- Moreover, the 
modular invariant partition functions may be classified, which corresponds 
in a certain sense to a classification of CFTs themselves. It has also been 
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shown that there is a deep correspondence with the ADE classification of 
simply laced finite dimensional Lie algebras [15|, 16]. See |10| and references 
therein for more details. 

1.4.2 Free boson 

It is now necessary to give a brief summary of the free boson theory which 
is a c = 1 CFT of a scalar field <j)(x, t) compactified on a circle of radius R. 
The Lagrangean is 

1 f L 

£=— dxd^(x,t)d^{x,t) (1.57) 

where L is the spatial volume, i.e. the theory is defined on a cylinder of 
circumference L. The field is quasi-periodic in the space direction, that 
is <p(x + L,t) = 4>(x, t) + 27rmR with m G Z. The integer m specifies a 
topological class of configurations obeying the above periodicity condition. 
Therefore one can characterize several sectors labelled by a pair (n, m), where 
-j| is the eigenvalue of the total momentum and m (winding number) is the 
eigenvalue of the topological charge, which are defined respectively by 

n = — / dxd t <p{x, t) , Q = — — / dxd x <f){x, t) 
Adopting the variables z and z, the field is expanded in modes as it follows 



0, z) = 0o - ip+ log z - ip- log z + % ^J - [ot k z k + a k z k ) 



k 



1.58) 



where (fio is the boson zero mode, p± — (^ ± |mi?j and 

[a n , a m ] = [a n , a m ] = nd n+mfl , [a n , a m ] = (1.59) 

The Virasoro generators take the form 

1 oo 1 oo 

L n = - £ «„-,«, , L n = - 5^ «„_ fc « fc (1.60) 

k=—oo fc=— oo 

The products above have to be considered as normal ordered. The modes a n 
and a n are annihilation operators for n > and creation operators for n < 0. 
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For any sector, the whole space of fields is obtained starting from the highest 
weight states \n,m) created from the vacuum by the vertex operators 

V {n , m) (z,z) = e^^ +p -^ 

of conformal dimensions Ar n ,m) = P+/2 an d A( nm ) = p 2 _/2, with successive 
application of creation operators on the states. Schematically one has 

ft = 0ftn,m<8>ftn,m (1-61) 

n,m 

where 7i„ im and 7Y„ jm are representations of Heisenberg algebra. Of course 
the states \n,m) are also Virasoro highest weights, therefore one could build 
up the whole Hilbert space by acting with the L n 's, but this would be more 
complicated. 

The boson Hamiltonian is expressed in terms of Virasoro operators as in 
Eq. ( 11.491 ) with central charge c = 1. 



1.5 Boundary conformal field theory 

We want to consider now a quantum field theory defined only in the half 
upper plane v > (see Fig. ( 11.21 )). One would like that the presence of the 
boundary on the axis u does not break several symmetries of the theory. In 
particular one expects that it is invariant under global rotations, dilatations 
and translations preserving the boundary and that these invariances should 
be implemented locally, i.e. conformal invariance realized in presence of 



boundary [17] 



A condition for boundary conformal invariance is that T 12 = T 2 \ on the 
axis u, which means that there is no flux of energy through the boundary. 
As a consequence the holomorphic and antiholomorphic components of the 
energy-momentum tensor are not independent, but T(z) = T(z) for 3mz = 0. 
So in the region 3mz < one can define the energy-momentum tensor 

T(z) = f(z) , 3mz < (1.62) 

Therefore in a boundary conformal field theory (BCFT), instead of having 
holomorphic and antiholomorphic parts of the fields in the half plane, it is 
possible to describe the whole theory with the holomorphic (antiholomorphic) 
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(u2,v2) O 



g ( ul ' vl > 



Figure 1.2: Semi-infinite geometry. The physically accessible half-plane 
is the upper one. 



part of the fields in the plane. 
Hamiltonian becomes 



For instance, in radial quantization, the 



H 



T(z)dz + c.c. 



1.63) 



c 



which means that the Hamiltonian generates the motion just outwards from 
the origin in the upper half plane. Because of the definition ( 11.621 ), the 
contour in Fig. ( 11.31 ) becomes a closed contour integral of T only: therefore 
the Hilbert space of a BCFT is described by a sum of representations of only 
one Virasoro algebra: 



W = 0«A) 



(1.64) 



Such a system is naturally represented mapping the half complex plane onto 
a strip of width L with some boundary conditions through the coordinates 
transformation w = (L/tt) ln(z). Note that, because of the single Virasoro 
algebra, the Hamiltonian in the cylinder now reads 



H 



z<* 



c 
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The space depends on the boundary conditions, it is organized into a vector 
space by a single Virasoro algebra strictly related to the boundary. Physically 
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Im z>0 




Im z<0 



Figure 1.3: Geometry of the integration contour for the boundary case. 



this implies that the same observables are associated to different representa- 
tions of Virasoro algebra in the bulk and boundary cases. 

1.5.1 Partition function and boundary states 

To study boundary conditions preserving conformal invariance and boundary 
states, it is useful to consider the BCFT on a cylinder within two equivalent 
quantization schemes, one in which the time flows around the cylinder, an- 
other one in which the time flows along that. In the first case the Hamiltonian 
depends on the boundary conditions a and b on the edges of the cylinder. In 
the second case, the boundary conditions are embodied in initial and final 
boundary states \a) and \b) and the Hamiltonian is obtained from the whole 
complex plane. It is convenient to introduce a partition function 



r exp ( — —H ab 



Tr q Hab 



1.66) 



with R and L respectively the circumference and the length of the cylin- 
der, q = e 2mT , t = iR/2L. In particular the cylinder can be obtained first 
mapping the upper half plane into an infinite strip, w = -log 2, and then 
imposing periodic boundary conditions in the R direction. The conformal 
invariance imposes that the spectrum of the Hamiltonian H^ falls into rep- 
resentations of the Virasoro algebra. If one calls M % ah the number of copies of 
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Figure 1.4: The same domain in different coordinates: the upper half 
complex plane punctured at the origin, the infinite strip with 
coordinate w = (L/2ir) log z, the circular annulus with coor- 
dinate u = exp(— 2iriw/R). 



representations in the spectrum, the partition function can be written 

Zab(q) = ^2KbXi(Q) 



as 



1.67) 



where the x«(qO' s are the Virasoro characters of the representation i. Inter- 
changing now the roles of L and R, that corresponds to the modular trans- 
formation r — > — 1/r, it is possible to regard the partition function as a trace 
of a Hamiltonian generating translation along the spatial coordinate. This 
can be done via the coordinate transformation u = exp(—27iiw/R), which is 
a plane distinct from the upper half z-plane (see Fig. ( jl.41 ) for a pictorial 
description of these transformations). 
The Hamiltonian reads now 



On the w-plane the boundary conditions are imposed propagating states from 
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the initial boundary state \a) to the final one \b). The partition function 
becomes 

Z ab {q) = (a\ e LR \b) = (a\ {f 2 ) ^^^ \b) (1.69) 

where q = e~ 2m ^ T . Under these coordinates transformations, the condition 
( 11.621 ) implies on the energy-momentum tensor 



u 2 T(u) = u 2 T(u) 



\u\ 



1, e 



2-kL/R 



This gives a strong condition on the boundary states \a) and \b) 
(L n - L. n ) \a) = , (L n - Z_ n ) |6) = . 



A solution to these equations is given by the Ishibashi states [18 

|j) = ^|j;iV)®|j;iV) . 



1.70) 



;i.7i) 



A ; 



where \j;N) and \j,N) are orthonormal basis of Virasoro representations j 
(N labels the different states within that module). Every boundary state 
should be a linear combination of vectors ( 11.711 ). For the partition function 
one has 



z ab(q) =^2(a\j)(j\b)Xj(q) 



1.72) 



Making use of the modular properties of the characters and comparing the 
expressions ( 11.671 ) and ( 11.721 ) it is also possible to give an exact valuation of 
numbers J\f* b . 

In the case of the free boson, the boundary states have to satisfy the 
stronger constraint (remember the form of Virasoro generators in the free 
boson scheme) 

(a n ± a n ) \a) = , (a n ± a n ) \b) = 

corresponding to Neumann and Dirichlet boundary conditions (T12 = 0). In 
particular the negative sign can be solved by states of type 



I a) oc exp 



E 

71=1 



Q.— n Q.— n 



n 



|0,fc) 
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therefore, any boundary state is a linear combination of vectors of this type. 
If one looks at this system on the cylinder and considers two boundary fields 
$ + and $_ acting on boundaries of both sides of the cylinder, one can calcu- 
late the partition function of the system explicitly (here we consider Dirichlet 
boundary conditions) and find 



Z 



DD =0C 



Z^' 



[oo 
n=l 



: ((J> + -(J)_+27rnr) 2 



;i.73) 



i 

r){q) 

where 4-y/vrr = R and 77(g) = g 1 / 24 \ 

For Neumann boundary conditions, or, equivalently, Dirichlet boundary 
conditions on the dual fields §±, 



q n ) is the Dedekind function. 



Z NN oc — — VV J 



$+ — &-+n/r) 



;i.74) 



These expressions have simple interpretations: one has to sum over all the 
sectors where the difference of hights between the two sides of the cylinder is 
$ + — $_ + 2irnr (Dirichlet) and $ + — $_ + n/r (Neumann). For each such 
sectors the partition function corresponds to the product of a basic partition 
function with hights equal on both sides, times the exponential of a classical 
action. The boundary states are, respectively, 



\B D ($ ± )) oc Y^ ex P 



ik$ 



± 



cxp 



-E 



L-t IT) *--*• — 



n=l 



??. 



|0,fc) 



B N $+ ) oc 



00 
y exp ( — 2ml<&±r ) exp 



i=— 00 



a-r,a. 



EVL- n Vi 
n 



n=l 



|Z,0) 



Details on this subject can be found in [19[ . 
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Chapter 2 



Integrable theories and finite 
size effects 



In this Chapter the main properties of integrable field theories are illustrated 
and the physical phenomena coming from finite size are introduced. 



2.1 Near the critical point 

It has been noticed that the analysis of the universality classes of two di- 
mensional statistical models includes the study of CFT living on the RG 
fixed points or, equivalently, from a statistical point of view, second order 
critical points and the description of the behaviour around them. The way 
a CFT behaves in a neighbourhood of RG fixed points can be observed con- 
sidering that any RG trajectory flowing away from such fixed points can be 
described, as briefly illustrated in the previous Chapter, by combinations of 



relevant fields $j, i.e. conformal fields characteristic of the given CFT [20 
The off-critical theory is described by the action 



A = A c ft + J2 Xi ^ x $ ^ x ) ( 2 -^ 



where Acft corresponds to the conformally invariant action and the Aj's are 
coupling constants. The Lorentz invariance in two dimensions is equivalent 
to rotational invariance once a Wick rotation is performed. Therefore, to 
have Lorentz invariance, one has to require the fields $, to be scalars under 
two-dimensional rotations, as a result their spin is s = A.j — \ = 0, which 
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implies that A, = Aj and the anomalous dimensions are 2Aj. The relevant 
operators $j in ( 12.11 ) do not affect the short distance behaviour of the theory 
because they are superrenormalizable, but they modify it at large distance 
scales. In particular, one can expect that through the RG flow the system 
either can reach another fixed point, in case described by another CFT, or go 
to a non critical point, hence it will correspond to a massive quantum field 
theory. In this scheme, the knowledge of the properties of the RG flow and 
of the theories linked by that is one of the main aims. For an unitary theory 
there exists a theorem, due to Zamolodchikov, which allows one to have an 
exact function describing the RG flow and, in the fixed points of such a flow, 
coinciding with the central charge of the underlying CFT. 

2.1.1 The c-theorem 



The c-theorem |21f] states that, given an unitary 1 + 1 dimensional quan- 
tum field theory endowed with rotational invariance and conservation of the 
energy-momentum tensor, there exists a function C(Aj) of the couplings Aj 
which is non-increasing along the RG flow and stationary only at the RG fixed 
points where it coincides with the central charge c of the corresponding CFT. 
One can demonstrate the theorem in a very simple way. Let T, and T be 
respectively the components with spin 2, and —2 of the energy-momentum 
tensor. The two-point functions among them can be written as 

, m/ , m , N1 F(mzz) 
(T(z,z)T(0,0)) 

(T(z,z)Q(0,0)) 

(e(*,z)e(o,o)> 



2 1 



G{mzz) 
z 3 z 



H(mzz) 

Z 2 Z 2 



where m is a mass scale parameter and F, G and H some scalar functions. As 
a result of the energy-momentum tensor conservation one deduces differential 
equations for the functions F, G and H 

F + - (G - 3G\ = , G-G+-(h-2h\=0 (2.2) 
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where F = dF(x)/dlogx. Defining the scalar function 



C = 2F-G--H (2.3) 



it follows 



C = --H (2.4) 

4 

The unitarity conditon of the quantum field theory imposes that if is a 
positive quantity, thus C is not increasing. At the critical point the energy- 
momentum tensor is traceless: G = 0. As a consequence, at fixed points, 
G = H = and F — | and hence the function C reduces to the central 
charge c. Moreover it is possible to extract a sum rule for the total change 
in the function C from short to long distances. In fact one can consider a 
CFT simply perturbed by one relevant field $ with coupling A and scaling 
dimension x = 2A. The trace of the energy-momentum tensor follows by 
direct calculation 

9(x) = 2ttA (2 - x) $(x) (2.5) 

and using the c-theorem, the total change of C is 

SC = 3ttA 2 (2-x) f dx 2 \x\ 2 ($(x)$(0)) (2.6) 

which means that one has an exact computation of the total change of the 
central charge along two different fixed points of the RG trajectory. If the 
CFT Ajj FT , from which the RG flow starts by application of a relevant scalar 
operator $, is known, in principle the CFT Aq FT in which the RG flow ends 
is also known. 



2.2 Integrable models 

An integrable quantum field theory (IQFT) is characterized by the existence 
of an infinite number of conservation laws, i.e. by an infinite set of con- 
served currents. In two dimensional systems one can use holomorphic and 
antiholomorphic indices and, in particular, as noticed before, an important 
characteristic of CFT in two dimensions is the decoupling of z and z depen- 
dence. So, in a CFT one can take as current any independent operator in 
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the conformal block of the energy-momentum tensor T(z, z) and, as a result 
of the conservation laws 

d g T(z,z) = d g T(z,z) = 0, (2.7) 

it follows that the CFT's possess an infinite set of conserved charges, i.e. 
they are integrable. It can be demonstrated that, given a polynomial T n (z) 
of the energy-momentum tensor T(z) at level n (and, equivalently, for T(z)) 
there are infinite integrals of motion of the form 

I n+1 = A dz T n (z) (2.8) 

which are non-trivial for n even. In the deformed theory defined in ( 12. 11 ), this 
set of conserved quantities - together with the decoupled variables z and z - 
is in general lost. Anyway, there exist the so-called integrable deformations 
of CFT in which an infinite number of conserved quantities survives and the 
theory can be treated non-perturbatively. 

2.2.1 Conservation laws 

For integrable theories originated by a perturbation of a CFT, the integrals 
of motion can be interpreted as deformations of the conformal conservation 
laws. For simplicity let us consider just one perturbing field $. 

In fact, if the equation d 2 T n (z) = for the linear combination T n (z) 
of fields in the conformal block of the energy-momentum tensor introduced 
above holds at criticality, this is not more valid off-criticality. The equation 
is deformed by several contributions 

d- z T n (z) = \Q^ X + A 2 QH + . . . (2.9) 

It is reasonably correct to think that the fields Q n \ have a limit in the 
unperturbed CFT, defined by some combination of fields in the Virasoro 
algebra of the CFT, with conformal dimensions (A 1 -", A^). The following 
relations hold 

AW = ra - Jfe(l - A) , A (fe) = 1 - fc(l - A) (2.10) 

with A the conformal dimension of the perturbing field $. This requires 
that the series ( 12.91 ) contains a finite number of fields Q n i\, as the conformal 
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dimensions of a CFT are bounded by below. In the simplest case there is only 
one field with conformal dimensions given in ( 12.101 ), hence one can consider 
that only the field Q n \ of first order in A appears, it has to be a secondary 
of the field $ with spin s = n — 1 and conformal dimensions (n — 1 + A, A). 
It can be demonstrated that 

d- z T n - d z e n . 2 = (2.11) 

where O n _2 is a field in the set of secondary descendant fields of $. This 
is an argument of the existence of perturbed conserved currents in IQFT. 
A simple example to show how to calculate exactly the conserved charges 
can be done just considering a minimal model 5DT r)S perturbed by a primary 
scalar field <& r>s (z, z). The action of the system is, following Eq. ( 12.11 ), 

A = A m + X J d 2 z <$> rjS (z, z) (2.12) 

The simplest conserved current which one can work with is the energy- 
momentum tensor T(z). Its OPE with the perturbing field is known from 
Eq. ( 11.331 ) and explicitly reads 

T(z, z)$ r , s {w, w) = — — + + . . . (2.13) 

[z — w) 2 z — w 

Considering now that the Ward identities can be written in terms of the 
conformal ones as 



(T(z,z)...) = (T{z,z)...) m + 

+ A f dwdw (T(z, z)® r , s (w, w)...) + (A 2 ) , (2.14) 

one easily arrives to the perturbed counterpart of the conformal conservation 
law at the first order in A 

d- z T{z, z) = A [(1 - &)d g $ r , a (z, z)\ (2.15) 

Therefore, identifying Q(z,z) = A(l — A)Q T)S (z, z), the first integral of motion 
is 

h= <bdzT(z,z)+ <bdzG(z,z) (2.16) 
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It is also possible to give several arguments to find higher conserved quantities 

In B. 



2.2.2 S-matrix approach 

In the analysis of a massive IQFT the scattering description has a very im- 
portant meaning. In general one can assume to have massive particles distin- 
guished by some label a, with mass M. a - Their momenta p% can be written 
in terms of a rapidity variable 9: E = _Mcosh#, P = .Msinh#. In any 
scattering process there are the "in-states", corresponding physically to a 
set of particles incoming from the infinite past, arranged by decreasing order 
of rapidities, formally described by \$i, . . . , On^n ' aN ■ At infinite time in 

a' a' 

the future, there are "out-states" \9[, . . . , #jv)o£/t n describing in general a 
different set of particles arranged by increasing order of different rapidities. 
The "in" and "out" states form a complete set of states of a local quantum 
field theory and they are connected by the operator called S matrix. 

In a IQFT, the existence of infinitely many conserved quantities has very 



important consequences [22] . In any scattering process it results that: 



• the number of particles is conserved, more precisely the number of 
particles of same mass is conserved; 

• the set of final momenta coincides with the set of initial momenta, 
{Pi 5 £>2? • • -Pn} = {p'iiP^i ■ ■ -P'n\i sucn that the scattering processes 
which take place are purely elastic. 

From this it follows that the S matrix factorizes into a product of 2-body 
scattering S matrices. To verify this, let Ai(9) be a set of non-commutative 
operators which describe the corresponding particles. They are regarded as 
generators of the infinite-dimensional algebra given by all possible products 
of type A ai (#i) A a2 (9 2 ) . . .A an (9 n ). The incoming and outgoing asymptotic 
states correspond to a decreasing and increasing arrangement of the rapidities 
9i. In this approach the S matrix is an operator such that 

Ai (9,) A, (9 2 ) = S* (0 12 ) A t (9 2 ) A k (9 t ) (2.17) 

where the relativistic invariance constrains the S matrix to depend on the 
difference of rapidities % = 9i—9j. These commutation relations are required 
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Figure 2.1: Factorization of the scattering. 



to be compatible with the algebraic associativity, which translates to the 
Yang-Baxter equation (YBE) 



okik2 
°i\i2 



7 12 



qjijs 
°k 1 k 3 



'13 



CJ2k3 
°k 2 i 3 



723, 



ok2k3 
12J3 



7 23 



qkiis 
°hk 3 



'13 



OJ1J2 

°k l i 2 



7 12J 



(2.18) 



It can be described by the following physical argument. One can consider just 
a particle, say a±, with its total momentum p±. Since the S matrix conserves 
momenta, one can think to conjugate it by the operator exp^p^x^), in fact 
it does commute with this operator, therefore this operation does not change 
anything. But it is a non-trivial physical action: conjugating the S matrix 
with exp(zp^x /i ) changes the space-time coordinates of the particle a\. In 
particular, choosing appropriately p±, one should be able to arrange for a\ to 
scatter with the other particles only when they are separated, therefore, the 
complete scattering process can occour as a succession of two particle scat- 
terings. So, proceeding inductively for any particle, the S matrix factorizes. 
Of course, for reasons of consistency, the scattering must be "associative", 
i.e. the scattering of three particles has to be decomposed into three pair- 
wise scattering with a result indipendent of which particular decomposition 
is used. From this the YBE follows and the Eq. ( 12.181 ) can be understood 
looking at the Fig. ( 12.11 ). 

The YBE is the fundamental constraint of the algebraic approach to 
IQFT's. The S matrix is not only characterized by the fact that it has to 
solve the YBE; there are several other requirements, as unitarity and crossing 
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symmetry 



Y, s T^) s lS-e) = ^ , sH(e) = s»(i*-$) (2.19) 



There exists also a consistent bootstrap principle, in the sense that there 
is a consistent set of couplings among the particles, signalled by poles in 
the S matrix at certain imaginary relative rapidities. These may be used to 
relate the S matrix elements to each other and, in a certain sense, to give a 
physical description of the particle system. In fact the S matrices are analytic 
functions in the complex plane of the Mandelstam variable s = (pi + p%) 
with branch cut singularities at {M.\ — M.2) an d {M.\ + M.2) > hence SM(9) 
are meromorphic functions of the rapidity. Usually the bound states are 
associated to simple poles, and in the bootstrap approach they are associated 
with some of the particles appearing in the asymptotic states. In particular, 
if 9 = iafj is a pole in the scattering process of particles a«, a,-, tha mass of the 
bound state a n is M.\ = M? + M.? + 2MiM.j cos afj. If the system has a non 
degenerate mass spectrum or presents a spectrum with degenerate particles 
which can be distinguished by their higher conserved charge eigenvalues, 
the S matrix is diagonal, therefore the YBE is satisfied trivially and the 
unitarity and crossing relations simply become Sij(9)Sij(—9) = 1, Sij{in — 
9) = Sij(9) where i denotes the antiparticle. In this case the bootstrap 
principle is encoded in the following functional equations 

Su(9) = Sv (9 + a* ) S lk (9 d - a\ k ) (2.20) 

where afj = 11 — afj. A pictorial representation of the bootstrap is illustrated 
in Fig. (QM . 

We conclude this section mentioning that there are also arguments to find 
the spin of the conserved quantities of an IQFT, by solving some particular 
equations encoding conserved charges eigenvalues, their spins and the angles 
ctij themselves with the help of special sets of resonance angles a^ . 

2.3 Integrable theory with boundaries 

At this point it is very useful to give a general idea of integrability when the 
two dimensional field theory is restricted to a half-line, or to a segment of the 
line (in space direction). It is clear that field theory defined in finite space 
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Out 




J k J 

Figure 2.2: Bootstrap relation 




with non trivial boundary conditions turns out to have more constraints than 
theories in infinite space or with (anti)-periodic boundary conditions. 

For example, suppose an integrable field theory describes a collection of 
distinguishable particles. It is not trivial to establish what is the spectrum 
of particle energies if the system is enclosed on the interval [—L,L]. Or, if 
the theory is confined to the space region x < 0, one might expect that the 
boundary at x = affects the particles approaching from the left. In fact, 
one has the almost minimal effect of reversing all the momentum-like con- 
served quantities and the preservation of all energy-like conserved quantities. 
It might be expected that the "out" state consisting of a single particle is 
proportional to the "in" state with the momentum reversed. 

The principal idea to solve the problem of integrability in the presence of 
the boundary (suppose, for the moment, it is at x = 0) is that particle states 
continue to be eigenstates of energy-like conserved charges and any initial 
state containing a single particle moving towards the boundary will evolve 
into a final state of a single particle moving away from the boundary [23 
i.e. 



\a,0) 



OUT 



K ab 16, 



'IN 



(2.21) 



where the states a and b correspond to multiplets of particles and K a b is a 
matrix which may mix the particles as a result of the reflection from the 
boundary. The first task in a boundary field theory is to determine the 
reflection matrices K for a given boundary condition. 

In what follows one assumes that the boundary does not affect particles 
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Figure 2.3: Boundary Yang Baxter equation. 



moving towards it until they are very "close" to it, therefore for all particles 
"far" from the boundary all the results of integrable theories hold. A first 
reasonable way to find out reflection matrices is to suppose that, given an "in" 
state with several particles, the order of the individual scatterings from each 
other and reflections from the boundary is irrelevant. If it is also supposed 
that these events take place factorizably, one obtains a boundary Yang Baxter 
equation (BYBE): 



K a (9 a )S ab (6 b + 6 a )K b (0 b )S ab 

= S ab (6 b -6 a )K b (6 b )S, 



<ib 



e n )K n 



(2.22) 



(see the Fig. ( J2.3J ) for a pictorial representation). 

Assuming now that the family of all masses and couplings remains the 
same in the presence of the boundary, the boostrap implies relations between 
the various reflection factors. Therefore, for example, if at special values of 
difference of rapidities ab the particles a and b form a bound state c, one 
might think of either a, b separately reflecting from the boundary in advance 
or after the bound state forms, or the particle c reflects from the boundary. 
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Figure 2.4: Reflection bootstrap equation. 



The reflection bootstrap reads 



K c (9 C ) = K a (9 a ) S ab (8 b + 9 a ) K b (9 b ) (2.23) 

where 9 a = 9 C - i9\ c and 9 b = 9 c + i9 a hc (see also Fig. (23}). 

One might also think to the possibility of bound states involving a particle 
and the boundary, i.e. the boundary can be excited [24|, [25], [26], this should 
be indicated by the presence of poles in K a {9 a ). For instance, a particle 
a and a boundary, that for simplicity here we denote with a, have, say, a 
reflection factor K^\ a pole at 9 a = iip^a can be interpreted as a boundary 
bound state (3. A two particle state a, b could also have the pole 9 a , in this 
case the particle b should be regarded as either reflecting from the boundary 
/3, or reflecting from the boundary state a, therefore it scatters twice with 
the particle a, once before and once after its reflection from the boundary 
state a. This implies a factorization assumption that algebraically is 



Kl (9 b ) = S ab (9 b + itfj K? (**) S ab (9 b - i< J (2.24) 

Pictorially, the boundary bound state bootstrap is represented as in Fig. 
( j2.5[ ). Finally one can also have unitarity and crossing relation as in the case 
without boundary 

K a (9)=K~\-9) , K a (9-i 7 ^jK- a (9 + i'^)S- aa (29) = l (2.25) 
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Figure 2.5: Boundary bound state bootstrap. The dashed line on the 
boundary represents pictorially an excitation of the boundary 
itself, i.e. the boundary bound state. 



In terms of lagrangean formulation of a quantum field theory with boundary, 
one has to consider that the first step is to take a Lagrangean evidently 
modified due to the presence of the boundary. 

For example, one can take for semplicity a field theory of a scalar field 
4> defined on the negative half line with a boundary at x — 0. If Co is 
the Lagrangean of the corresponding IQFT without boundary, the modified 
Lagrangean might take the form 



Lb = Step(— x)£q — 5(x)B((p) 



(2.26) 



Therefore, the field equations are restricted to the region x < with a 
boundary condition at x = 



0, 



dB{(f>) 



(2.27) 



The main problem is to establish what value of B is compatible with inte- 
grability. In particular, one has to take into account that the presence of 



44 



2.4 -- From Integrable models to critical points 



the boundary removes the translational invariance, it is not expected that 
the momentum should be conserved. The energy should be. One can expect 
that energy-like charges might continue to be conserved, while momentum- 
like charges cannot be. The way to find integrability in presence of boundary 
is, essentially, to generalize the standard Lax pairs approach including the 
boundary. Moreover, it is also requested, in order to have a complete de- 
scription of a boundary IQFT, an inverse scattering procedure. It does not 
exist a systematic determination of integrability in presence of boundary and 
the full set of reflection factors for any specific integrable model is in gen- 
eral unknown. Anyway, many progresses have been done in the search for 
integrability and knowledge of reflection factors in the boundary integrable 
theories as sine-Gordon and affine Toda theories 127, 28, 29 . 



2.4 From Integrable models to critical points 

A link between IQFT coming from perturbations of a given CFT and fac- 
torizable scattering theories can be done by studying the finite size effects 
associated to theories due to the finiteness of space regions. What one usu- 
ally works with is the scaling properties of a quantum field theory and, in 
particular, the scaling functions and their behaviour along the RG flow. For 
example, it can be useful to deal with special quantities, as the energy of a 
system, whose scaling in particular space regions of definition is well known. 
It has been noticed that, in a cylinder, a CFT develops a sort of Casimir 
effect due to the finiteness of the spatial extension L, therefore, using for- 
mulae ( 11.481 ) and ( |1.49[ ), it follows that for a conformal state \i) the energy 
eigenvalues behave as 

Ei = ~~d with c * = c - i2 ( A *+ A ( 2 - 28 ) 

Of course, off-criticality this dependence of the energy on L might change. 
However, all the variations of the dependence can be put in a dimension- 
less function such that it reproduces ( 12.281 ) in the limit L —> because the 
dependence 1/L is fixed by dimensional reasons. This allows to introduce 
the dimensionless parameter / = M.L and the scaling function q(/) which 
describe the behaviour of the energy state \i) in the regions I — > and I — > oo 
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through a sort of generalized version of the Eq. ( 12.281 ) 

El (l) = -^ (2.29) 

with the request that q(/) — > c — 12 (A, + Aj) for I — > 0. The presence of 
a mass parameter .M has not to be surprising because the function q(7) is 
describing the system flowing away from the critical point, it is thus natural 
to introduce a mass scale. Moreover, in this scheme, the IR and UV limits 
take a stronger physical meaning as well. The limit / — > (UV) can be also 
interpreted, taking the spatial extension L fixed, as a limit for very small 
masses, i.e. a very large energy scale system in which the specific masses 
have no more significance. This situation, in fact, is usually referred as UV 
(high energies) limit in quantum field theory. On the other hand, the limit 
/ — ► oo can be seen as a physical limit where the masses are very large, 
acquiring, thus, a great relevance, i.e. one realizes a theory with several 
scattering massive particles. Let us notice that for a BCFT the relation 
(07281) becomes 

Et = with Ci = c- 24A; . (2.30) 

24L K J 

One way to compute the scaling functions is to make use of perturbation 
theory, and in our case of conformal perturbation theory. However, it is 
valid only in a neighbourood of a fixed point and, only sometimes, it allows 
to make contact between different fixed points, i.e. CFT's. The real aim is 
to link CFT to massive scattering quantum field theory, which are evidently 
not conformally invariant. This leads to the application of non-perturbative 
methods. These methods are essentially: 



the Truncated Conformal Space Approach (TCSA) [30] consisting in di- 



agonalizing the truncated Hamiltonian numerically. Even if it is non- 
perturbative and applicable also to non-intergrable systems, it is af- 
fected by numerical errors and it does not give any analytic control of 
the scaling functions; 



the Thermodynamic Bethe Ansatz (TBA) |3l| based on the idea of 
implementing the thermodynamics of a statistical system of particles 
interacting via a given S matrix. In particular, the TBA allows to cal- 
culate exactly the free energy of the system and other thermodynamical 
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quantities as the entropy. Once the role of space and time have been 
exchanged, the finite size effects of the vacuum energy and first excited 



states 



32 



33] energy appear and can be studied analytically; 



the Nonlinear Integral Equation (NLIE) [34, [35], 36], |37j, which is an 
approach allowing to obtain scaling functions for a given IQFT starting 
from its definition as a quantum field theory regularized on a finite 
lattice. 



In the next Chapter we will be introducing this last method for a particular 
integrable system and giving a description of its finite size scaling between 
the UV and IR limits. 



47 



Integrable theories and finite size effects 



48 



Chapter 3 



Boundary Sine-Gordon Theory 



The sine-Gordon model in a strip with Dirichlet boundary conditions is in- 
troduced. Particular emphasis is given to its regularization on the lattice via 
an XXZ spin chain with magnetic fields at the boundaries. Bethe Ansatz 
equations and energy eigenvalues are derived as well. 

3.1 General results 

Quantum field theory with boundaries has been for many years an interest- 
ing subject of theoretical investigations and has been found to have many 
applications too, for example, in dissipative quantum systems [38], [39J. In 
particular, boundary IQFTs provide some useful tools as to study the struc- 
ture of the space of boundary interactions, which have great significance in 
open String Theory |40|, as to probe some fascinating physical phenomena in 
Condensed Matter Physics, for instance, the Kondo effect and the quantum 



Hall effect M, \42 



Let us consider the sine-Gordon model 

A SG = / dt / dxa(<f),d^) (3.1) 



+00 p+oc 








dt / 


dx a 


(0, 


d,cj>) 


00 J — 00 








= \ (<9,0) 2 


2 
m 


cos {(3(f)) 



with 

a (0, d^) 
where 4>(x,t) is a scalar field and (3 a dimensionless coupling constant. The 

49 



Boundary Sine-Gordon Theory 



field theory in the presence of a boundary, say at x = 0, can be defined by 
the following action 



/+oo r0 p+oo 

d£ / dx a (0, <9 M 0) + / dtB(</>(0,t),dt<f>(0,t)) (3.2) 

-oo J— oo J— oo 

where the possibility of boundary degrees of freedom other than the boundary 
value of the bulk field 0(x, t) has been discarded. It has been found in |24f 
that the boundary action preserving the integrability is of the form ( 13.21 ) with 

B (0(0, t), d t <f>{0, t)) = -ix cos (j (0(0, t) - O )J (3.3) 

where /i and 0o are f ree parameters. The boundary reflection amplitude has 



been derived by using boundary Yang-Baxter equation [24, I43J- Exept for 
the case /i — oo, the boundary value <p(x = 0, t) is not fixed in the boundary 
theory ( ]3.2l ) and hence the topological charge 

Q = ^J **^<l>(x,t) (3.4) 

is not conserved. We will refer to the theory ( 13.21 ) with B (0(0, t), 9^0(0, £)) 
given in ( 13.31 ) as boundary sine- Gordon model. 

Let us give a brief summary of the theory in the bulk, further details can 
be found in [22). The bulk sine-Gordon model is known to be integrable at 
both the classical and quantum levels [44j, 45]. The theory has an infinite 
number of degenerate vacua with the discrete symmetry — > + ^-m, with 
m G Z. The spectrum consists of a soliton a and an antisoliton a and a 
number of soliton - antisoliton bound states b n , n — 1, 2, . . . < -, usually 
called breathers with 



P=3T 32 ( 3 - 5 ) 



87T-/3 2 

The soliton and the antisoliton have the same mass A4 S , while the mass of the 
n th breather is A4 n = 2.M s sin (^f 2 ). The topological charge of the soliton 
(antisoliton) is 1 (—1), while the breathers are neutral. It is always conserved. 
As a consequence of the integrability, there is no particle production, i.e. the 
scattering is factorized and, if one denotes the soliton S matrix as S^(8) 
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with a, b, c, d taking the value + (— ) if the particle is a soliton (antisoliton), 
the non-zero amplitudes are S++(9) = SZZ(9) = a(9), which represent ei- 
ther the soliton - soliton scattering or the antisoliton - antisoliton scattering, 
S~\z{9) = Si +(9) = b(9), which denotes the soliton - antisoliton transmission 
and finally the soliton - antisoliton reflection Sli(9) = Sll(9) = c(9), where 
the functions a(9), b(9) and c{9) explicitly read 

r a(0)=sin(W)p(u) 

< b{9) = sin (A p{u) (3.6) 

c{9) = sin (j) p(u) 

with u = —%9 and p(u) can be written either in terms of Gamma functions 

as 



p(u) 



where 
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or in terms of the Barnes' diperiodic sine function £2 (x \ u>i,u>2) [46|, |47| , 
which are meromorphic functions parametrized by the pair of quasiperiods 
(u>i,u>2), as 



p(u) 



sm 



1\ 



P 



5*2 (tt — u I p7r, 27r) S 2 [u \ pn, 2tt) 

S'2(7r + ^ I J>7T, 27r) £2 (— U I p7T, 27r) 



(3. 



Let us mention that, due to the properties of Gamma functions, the function 
p(u) and, therefore, Eqs. ( 13.61 ) can be expressed in integral form. As an 
example, we write the soliton - soliton (antisoliton - antisoliton) scattering 
amplitude as 



a{9) = exp 



dk 



Ak9 



sinh |(p — \)k 



2%k 2 sinh ^pk cosh | k 



(3.9) 



51 



Boundary Sine-Gordon Theory 



this expression will be acquiring a great significance in next Chapters. The 
amplitudes b(9) and c{9) have simple poles at 9 = i(ir — mrp) and a(9) and 
b{9) have poles at 9 = inirp. They correpond to the creation of breathers b n 
in the forward and cross channel respectively. 

The model we are going to deal with is the special case with two bound- 
aries, on the left and on the right of a semi-infinite strip, infinite in time 
direction and from to L in space direction. The boundary conditions are a 
special restriction of the most general one ( 13.31 ) preserving the integrability, 
i.e. they are of Dirichlet types. We define the Dirichlet sine-Gordon (DSG) 
model in a strip of lenght L through the Action 



•Adsg = ~ / dt dx 

* J-oo JO 



nir 



(a i 0) 2 -(^) 2 + -|cos ; 
P 



(3.10) 



with the constraints fixing the value of the field at the boundaries: 



27T 2-7T 

0(0, f) = (p^ + — m_ and <j>{L, t) = <p + + — m + (3.11) 

P P 

with m± G Z. In the Dirichlet diagonal case the topological charge can be 
written as 

Q = i- J dx^ix, t) = -S- (0 + - </,_) + m + -m. (3.12) 

The model enjoies the discrete symmetry of the field (p — » (j) + %m and simul- 
taneously (j>± — > 4>± + ^jin. The charge conjugation symmetry is also guar- 
anteed provided — > —<p and (fi± -^ —(fi± simultaneously. It sends Q — > —Q, 
therefore, one can restrict attention to the study of positive Q and then 
act with this tranformation to obtain states with negative Q. Then, only 
processes conserving the topological charge are allowed, in particular, any 
particle of the theory scattering off the boundaries has to conserve its topo- 
logical charge. Observe that the periodicity allows to restrict the boundary 
parameters to the range < <f>± < ^. Scattering off the boundary processes 
are described by the diagonal term of the most general boundary reflection 
amplitude found in |24|. Therefore, if one introduces the reflection factors 
for the solitons V±(u) (a soliton or an antisoliton incident on the boundary 
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is reflected back unchanged) and Q±(u) (a soliton is reflected back as an 
antisoliton, or vice versa), one has that they take the values 

Q±{u) = (3.13) 



V±(u)=R (u 
where 



&(? TpS± + n-u \ P 7r,27r)S 2 (^TPt± + u \pk,2tt) 

(3.14) 



N g 2 (j- M |f,27r)g 2 (f + M |a 2 7T) 

° ^ 2 (f + n |ffi 2 7r)5 2 (f -« |f,27r) [60) 

£± = —n- (3-16) 

and S 2 (as | w\,w<i) are the Barnes' diperiodic functions introduced above. It 
has been shown |48| that there exists a quite intricate set of poles of the 
diagonal reflection factors V±(u) on the imaginary axis, which correspond to 
boundary bound states. In particular, it has been stated that the boundary 
bound states exist only for particular values of the topological charge and 
they are represented by the set of poles 



where 



{I'm, Wrm ■ X > Un i > Wm i > V n % > • • • > j (3.17) 

v n = p£,±- ^(2n+l)p with n = 0,1,2, ... (3.18) 

u! m = 7T — pi± — — {2m — \)p with m — 1,2, ... (3.19) 

In order to investigate all the RG flow of this theory, from the IR limit, 
i.e. the scattering theory limit, to the UV limit, i.e. the conformal limit, 
and associate states of the scattering particles interpretation to states of 
the c = 1 conformal field theory, we will construct the DSG model as the 
continuum limit of a lattice model by an inhomogeneus XXZ spin chain with 
fixed boundary magnetic fields, i.e. diagonal reflection amplitudes on the 
boundaries, in a Hamiltonian formalism. 
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3.2 Boundary XXZ model 

The homogeneous antiferromagnetic XXZ spin | model in a chain of N sites 
spaced by a, coupled to magnetic fields h- and h + at the left and right 
boundaries respectively, has Hamiltonian 



aH (h+, h-) = -jJ2 {«+! + ay n a l+i + cos 7«+i) + h-al + h+a] 



n=l 



(3.20) 
Here ct", a = x,y,z are Pauli matrices acting on the n-th site, so that 

< 7 < 7r, so that the isotropic case XXX is reproduced for 7 = 0. The 
parameter 7 is often referred as anisotropy of the chain. The model ( 13.201 ) is 
invariant under global U(l) rotations of the spins (a^ = ~ (a* ± a^)) 

—± > ^i* w ^-i z z 

°n -> e a n » ^n ~* °n 



conserving the third component of the total spin 

JV 



1 - 7 



n=l 



Therefore, states can be classified according to their value of 5 1 , which in a 



V JV 1 -1 AT 1 JV 

.e an values o = - 

The transformation 



chain of N sites can take all values S = — -j, — y + 1, ••-, y — 1, 2 



is another important symmetry, under which 

C- 1 H(h + ,h_)C = 'H(-h + ,-h_) and (T^C = -S 

which means that states of a Hamiltonian with given h + and /i_ with negative 
spin are equivalent to states with positive S of the corresponding Hamiltonian 
at the same 7 but with reversed boundary magnetic fields. In this sense one 
can decide to restrict to consider S > states only, as the negative ones will 
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Figure 3.1: Pictorial representation of the XXZ spin chain. The system 
is invariant under the trans fotmation C~ l TL{h + ,h^)C and 

c- l sc. 



be obtained by just changing sign to both the boundary magnetic fields in 
all formulae. 

However, unlike the case of a single boundary and a semi-infinite chain, 
this symmetry is not enough to allow to restrict to both h± > 0. Indeed it 
is physically evident that the case with parallel magnetic fields (h± > 0) is 
different from the antiparallel one (h + > , /i„ < for example). The only 
allowed inversion of magnetic fields is the simultaneous one. This symmetry 
will be related to the charge conjugation of the DSG model. In what follows 
we will consider the antiparallel case and choose /i_ > and h + < 0. 

3.2.1 Double row transfer matrix 



The Hamiltonian ( 13.201 ), as well as its inhomogeneous generalizations, can 
be constructed in a transfer matrix framework. Together with the two di- 
mensional spaces V n on which the cr"'s act, consider an auxiliary space Vo, 
isomorphic to the V n 's. Introduce the sl(2) q .R-matrix 



Rn 



1 + c(u) 



[u 



*-n*-m 



1 - c(u, 



a z n a z m + b(u) [cr+cT m + o n cr+] 
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and the quantum Lax operator 

L n (u) = P nfi R n ,o(u) 

where P n , m is the permutation matrix P„ jm = cr+<j m + cr~cr+. The w n 's are 
fixed parameters called inhomogeneities and 

, . . sinh u i sin 7 

b{u) = . — ■ , c[u) - 



sinh(u + i'-f) ' sinh(w + ij) 

Rn,m{ u ) satisfies the Yang-Baxter equation. 
Define the monodromy matrix 



JV-l JV-l 



T(«| U ) = n l m = n r ^-^ u - u ^ = ( cm dm ) 

which is a 2 x 2 matrix in the Vo space, whose elements are operators in the 
Hilbert space V = (S) n =i Hi> an d u — ( M i; •••? u n)- The trace in the Vo space 
of T 

t(u|u) = Tr T(w|u) = A(u) + B{u) 

is well known to be an operator in V called row to row transfer matrix. As 
a consequence of the Yang-Baxter equation, on a periodic chain 

[t(«|u),t(v|u)] = ,VU,!)GC 

In the case of non trivial boundaries, however, this row to row transfer matrix 
is not enough to define the integrals of motion. One has to introduce a more 
complicated object that takes boundary effects into account. Introduce also 
the boundary K- matrices K±(u) = K (u, h±) in Vq where 



K(u, h) = (cosh u) 1q + i (h sinh u) er ( 







which satisfy the boundary Yang Baxter equations. Define the 2x2 matrix 
in Vo 



A(u) B{u) 



U(u|u) = T{u\n)K+{u)T-\-u\n, ■ 

L[u) V[u 

In terms of this object the double row transfer matrix can be introduced 

56 



3.2 - - Boundary XXZ model 



d(u|u) = Tr K-(u)U(u\u) 

having the fundamental integrability property that [d(«|u), d(u|u)] = 0, 
Vu, v G C. In the homogeneous case u = 0, the Hamiltonian ( 13.201 ) is well 
known to be reproduced by: 



aH 



du 



logd(/u, 0) 



N 



u=0 



^L n (0) + -K + (0) 



71=1 



Tr K (0)L N (0) 



where dots denote derivatives with respect to u, and the label refers to the 



auxiliary space (for details see [49 



The bare continuum limit N — > oo, a — > while Na = 
is known to give the Hamiltonian of a massless free boson 



L remains fixed, 
i)(a;) compactified 



on a circle of radius R 



2(71-7) 
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Among the many possible deformations of the Hamiltonian ( 13.201 ) leading 
to sine-Gordon in the bare continuum limit, we choose the one introducing 
alternating inhomogeneities u n = (— l) n A in the sites of the chain. This 
construction has the advantage to preserve integrability of the lattice model, 
thus allowing to use the Bethe Ansatz techniques needed for our purposes. 
It has been known for sometime [51, 52] to give a correct construction of 
sine-Gordon model in the bulk in cylindrical geometry, when the appropriate 
scaling limit is chosen, with periodic or twisted boundary conditions. It is 
then natural to expect that the same construction in the presence of boundary 
magnetic fields h± can also provide an effective tool to define the renormalized 
DSG theory in a strip of lenght L with fixed conditions at the boundaries. 
It is worthwhile to recall that, as the homogeneous XXZ chain is equivalent 
to a 6-vertex model on a square lattice, this modified XXZ chain is also 
equivalent to a 6-vertex model, but - as a consequence of the introduced 
inhomogeneities - defined on a lattice rotated by 45°, i.e. on what can be 
thought as a Minkovski space discretized along the light-cone directions. This 
is why this construction is often referred as light cone lattice construction of 



the sine-Gordon model 51 



The double row transfer matrix, in this case [49], becomes 



d(u |u) = d (u, A) = K N (u, h + )T (u, A) K x (u, h_) f (u, A) (3.21) 
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where 



K N (u, h+) = Tr [R N>0 {u + (-l)^A) K{u + i 7 , h + )R N , (u - (-if A)} 



N-l 



T (U, A) = ] [ i^-n+l,iV-n (« " (-1)^"^) 



71=1 

AT-1 



7>,A) = n^n,™+i(w+(-l) n A) 

n=l 

It is possible to demonstrate that the transfer matrix ( 13.21J ) is an even func- 
tion of A up to an irrelevant multiplicative function of A only, indeed, sim- 
plifying the quantity K N (u, h+), one easily gets 

K N (u, h+) = f («, A) K N (u + i7, h+) (3.22) 

where / (u,A) = 1 + c [u — (— l^A) c(u+ (— l^A) is evidently a numerical 
function. Then define 

r (u, A) = K N (u + i7, fr+) T (u, A) i^ («, fc_) f (u, A) (3.23) 

and observe that d (w, A) = / (u, A) r (u, A). The following identity holds 

r- 1 (-u, A) r' (-«, A) = r" 1 («, A) r' (u, A) (3.24) 

In fact, given r (w, A) and making use of the fact that the matrix Kn (u + ij, h. 
is diagonal, one has 



r (u, A) = K N (u + i7, /i+) T (u, A) K x (u, h-) T (u, A) 

= T (u, A) K x (w, /i_) f («, A) fGv (« + 17, fc+) (3.25) 

Now, due to the properties of R matrix and K matrices 

Ri^-u) = Rr}(u) [K, (-u, h±) = iq 1 (u, h±) 

Ri,j(u) = Rj,i(u) [Kj (u, h±) = Kj l (-it - 2z'7, h±) 
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one performs the following computation 



r (-u, A) = K N (-u + ij, h+) T (-u, A) K x (-u, hJ) f (-u, A) 
= T (-u, A) K x (-u, h-) f (-u, A) K N (-u + 17, /i+) 
= [R N - ltN (-u - (-l^A) . . . R lj2 (-u + A) K x (-u, h, 
x [ Rl2 (_ u _ A) . . . ifo_ 1)2iV (-u + (-1)"-^) AT* (-« + 17, h+j\ 
= [R~ N \ N (u + (-l^A) . . . Rll (u - A) K^ (u, h. 



x 



X 



[^1,2 (« + A) . . . i^ 1 .^ (« 



\JV-1 



A) A^ 1 (it - 17, /i H 



[i#_ liJV (u + (-l^A) . . . R~l (u - A) K- 1 (u, h.)] x 



[#7J (u + A) ... R~\ N (u - (-lf~ l A) K N l (u + i 7 , k+)] 



Ktv (« + 17, M T («, A) A'i («, h-) f (u, A) 



and, taking the logarithm, 



log r (—u, A) = — log r (u, A) 



(3.26) 



which proofs the statement ( 13.241 ). This means that the logarithmic deriva- 
tive of the operator r (u, A) is an even function. So we have that 

|-logd(u,A) = A-hg f(u, A) + ^-\ogr(u, A) (3.27) 

because / (u, A) is a numerical function. At this point all the commuting 
integrals of motion can be defined in terms of d (u, A) as 



Q 



(k) 



a du k 



logd (it, A) 



(3.28) 



and, in particular, the Hamiltonian turns out to be 



i d 



H = *' = s" sdW) 



(3.29) 



u=A 



An explicit expression for our case with diagonal reflection matrices is 
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jv-i 



-iaH = J2 R n!n+1 (2A) Rn,n+1 (2A) + 

n=l 
N-l 

+ J2 ^i + i(2A)i?n-i,n(0) J R n , n+ i(2A) 



JV-1 



+ E ^, 1 „. + l( 2A )^+l,n+2(0) J Rn,n + l(2A) 



n=odd 



+ KZ 1 (A) K- (A) + K; 1 (A + vy) K + (A + vy) (3.30) 

For A — > the Hamiltonian ( 13.201 ) is reproduced. 

When A — > oo, N — > oo and a — > 0, while L = Na stays fixed, contact can 
be made, along lines similar to those in |53|, with the lagrangean formulation 
of DSG model. The XXZ anysotopy 7 is related to the SG coupling j3 by 
(3 2 = 8(7r — 7). It is often convenient to use the parameter p defined in ( 13.51 ): 
p= f -1, (0<p<oo). 
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3.3 Bethe Ansatz equations 

The Bethe Ansatz equations for the boundary XXZ chain ( 13.201 ) have been 
written some years ago |49|, 54J, using an algebraic approach. It is straight- 
forward to generalize them for Eq. ( 13.301 ) with the introduction of the al- 
ternating inhomogeneities. Eigenstates of the transfer matrix d(u, A) can be 
constructed applying repeatedly B(uj) operators to the ferromagnetic vac- 
uum |0) 

B(m) . . . B(u M )\Q) (3.31) 

The all dinstinct numbers U\,. . . ,um are called roots. They are in number 
of M, (M < N) and must satisfy the Bethe Ansatz equations 



[si/2 (#j + e) 81/2 & - e)] N s H+/2 (^-) s H _, 2 (-&J 



1 TV 
M 

H s 1 (# j --& k )s 1 (0 j +0 k ) (3.32) 
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where 



7T 



„ N sinh -($ + ivir) n 

s v (#)= . .*). — £ and = ~u,e=:~A 3.33 

smh - [v — ivk) 7 7 

The parameter H± is defined such that h± = sin 7 cot | (H± + 1) and we 
choose as fundamental domain — p — 1 < H± < p + 1. Notice that for 
H± = 0, i.e. h± = 1 + cos 7 = /i c , the boundary terms in the Bethe equations 



disappear and the theory becomes SL q (2) invariant [55, 56 



The antiferromagnetic vacuum turns out to be a set M ■ - y of real 
roots and it exists for N even only. In the region < 7 < n and for 
small enough boundary magnetic fields, this is the true ground state of the 
model. For N odd instead the states with lowest possible total spin have 
M = ^rp- roots and one hole. However, to deal correctly with the continuum 
limit one has to consider both N even and odd sectors. The symmetry of 
{$j} —> {— $j}, evident from the Bethe equations, implies that only roots 
with positive real part are independent parameters characterizing the Bethe 
states. The value dj = is a solution of Eqs. ( 13.321 ) for any N and M. 
However, the corresponding Bethe state would vanish, one has to subtract 
this unwanted root, i.e. to create a hole at 1? = 0. 

The domain of roots distribution can be considered as a semistrip U + of 
the complex $-plane 

( 22 

U + = <•& e C : 9te(0) > , -— < 3m(0) < — or 

I 2 7 2 7 

7T 2 



9Wtf) = , < 3m(tf) < — } (3.34) 

2 7 J 

2 

This excludes another unwanted root at i|- as well and considers only half of 
the imaginary axis, as it should for symmetry. For computational purposes, 
it is useful to double this strip by mirroring all the roots 

U = {■& e C : 9te(tf) e R , < 3m(ti) < — > (3.35) 

I 27 27 J 

where to each root dj is associated its mirror root $_j = — dj. Define the 
function 

<^(tf)=7r + zlog Si ,(tf) (3.36) 
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with the oddity condition tp „(—•&) = —<p v V&) fixing the fundamental branch 

2 

of the logarithm. It is periodic in the imaginary axis, with period i— and 
real on the real axis. We choose as fundamental peridicity the strip v G U. 
Singularities of this function appear along the imaginary axis: 

9te(tf) = , 2fai(0) = ±n(k(p + l)-v) , keZ (3.37) 

The fundamental analyticity strip is limited to |3ra('*9) | < 7rmin(z/, p + 1 — v). 
In terms of the function ( 13.361 ) the logarithm of the Bethe Ansatz equations 
can be expressed as 



N [ipx/2 (^ + 0) + <p 1/2 {dj - 0)] + <PH + /2 i'&j) + VH-/2 (^) + 

M 



k=-M 

3.3.1 Energy 

An useful result can be achieved from the diagonalization of the double row 
transfer matrix ( 13.211 ): eigenvalues of the Hamiltonian can be expressed in 
terms of solutions of Bethe equations ( 13.32) , due to the formula ( 13.29) . There- 
fore, given a Bethe state characterized by M roots in the set U, the corre- 
sponding energy turns out to be 



1 M d Id 

E = ~Y a £ de ^ 1/2 (e " ^ + Vl/2 (e + ^ + a de^ 1/2(0) (3 ' 39) 

fc=l 
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Nonlinear Integral Equation 



In this Chapter the fundamental Nonlinear Integral Equation for the Bethe 
Ansatz is derived following the approach introduced in [34J, [35], [36], 1371 ] . In 
the last section a general formula for the energy of the model is given, it 
allows to define and study scaling functions in the next Chapter. 

4.1 Counting function 

Define for iJeU the counting function 

Z N {d) = N [p l/2 (tf + 6) + <p 1/2 (0 - 6)] + <p H+ /2(#) 

M 

k=-M 

where the function (f^fi) has been defined in the last Chapter and the values 
i? fe) k = -M,-M + 1, . . M - 1,M, are solutions of the Bethe Ansatz 
equations ( 13.321 ). It is evident that in terms of the counting function Zjv($) 
the logarithm of the Bethe equations simply becomes 

Z N (fy) = 2-irIj. (4.2) 

The last two terms in Eq. ( J4.1J ) take care of the fact that in the second mem- 
ber of ( J3.321 ) it is not included the factor with k — j and that the root $o = 
is an unwanted solution and it has to be subtracted. Eq. ( ]4.2l ) plays the role 
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of a quantization condition for Bethe Ansatz equations roots, therefore, in 
the following, we will refer to the integers Ij as quantum numbers. 



4.1.1 Bethe roots 

The analytic structure of Z N ($) suggests to classify the roots and related 
objects of the Bethe equations as follows 



• real roots Xk, k = 1,2, ... Mr, Xk > 0: strictly positive real solutions 
of Eqs. Q332} and M ); 

• holes h-k, k — 1, 2, ... N H , h k > 0: strictly positive real soltutions of 
( M that are not solutions of ( 1X321 ); 

• close roots c k , k = 1,2, ... Mq- complex solutions of Eqs. ( 13.321 ) and 
( 14.21 ) with SHc(cjfc) > and imaginary part in the range < |3m(cfc)| < 
min(7,7r-7); 

• wide roots Wk, k — 1,2, ... Mw'- complex conjugate solutions with 
^(wk) > and min(7, ir — 7) < \3m(w k ) | < f • 

It is also useful to define M-jy, the number of wide roots with positive imag- 
inary part and M\^ those with negative imaginary part. Clearly, it holds 
Myy + Myy = M w , while M^y — M^y = M S c is the number of self-conjugate 
roots with 3m(wk) = |, whose complex conjugate is the root itself, due to 
the iir periodicity of Bethe Ansatz equations. 

The function Zjv($) for $ G K is globally monotonically increasing. How- 
ever, there may be points where locally Zn(i9) < 0. In particular, one can 
call special roots or holes s k the objects such that Z N (sk) < 0. If such a spe- 
cial object appears, there must be two other objects (either roots or holes) 
with the same quantum number as well. This is imposed by the global in- 
creasing monotonicity of the counting function. Moreover, as two roots with 
the same quantum number are not allowed in the Bethe equations, only one 
of this triple can be a root, the others are forced to be holes. The number of 
special objects will be indicated as N$. 
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w. 



Figure 4.1: Bethe roots and holes in the physical strip. The dashed 

2 

lines denote the periodicity — and the dotted lines the value 
±7rmin(l,p). 



4.1.2 Counting equation 

One may relate the numbers of various types of roots to the total third 
component of the spin of the system 



8-1-U. 



(4.3) 



In order to do that, one makes use of the asymptotic values of the functions 



lim <p u {$) = ±7r(7r — 2v) 



(4.4) 



which allows one to express the ±oo asymptotics of the counting function as 



Z N (+oo) =2tt (N-M) + (AM - 2N) 7 

- 2 1 H + 4(vr - 7) + 2vrSign(p - l)M^ 



(4.5) 
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Zjv(-oo) = - 2vr (N - M) + (AM - 2N) 7 

+ 2 1 H + 4(tt - 7) + 2vrSign(p - ljAfJ, (4.6) 

where H = + ~ . Substituting Eq. ( 14.31 ) in Eqs. ( 14.51 ) and ( 14.61 ), one gets 

Z N (+oo) = 2n(2S + M) + 2vra + 2vrSign(p - 1)M^ (4.7) 

Zjv(-oo) = -27r(25 + M) - 2na + 2vrSign(p - l)M ] w (4.8) 

with a = ^~ xj. Define 7 max to be the largest integer smaller than 

Zn(+oo)/2ti, and 7 m i n the smallest integer greater than Zn(—oo)/2ti. Evi- 
dently it holds 

/max = 2S + M+ [a\ + Sign(p - 1)M^ 



J min = -2S -M- [a\ - Sign(p - 1)M, 



T 



where \_a\ denotes the integer part of a. Therefore 

Zjv(+oo) = 27r/ rnax + 2tt (a - [«J ) (4.9) 

Z N (-oo) = 27r/ min - 2vr (a - |aj ) (4.10) 

Now, considering the total number of real Bethe solutions on the real axis, 
one has that 

/max - Imin = 2 (M R + N H - 2N S ) (4.11) 

and then, combining Eqs. ( 14.91 ), ( 14.101 ) and ( 14.111 ) and substituting the values 

°I -<max an d Jmin) 



N H - 2N S = 2S + M C + Step(p - l)M w + 

where 
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2 (p - S) H 



p+ 1 p+ 1 



(4.12) 



4.2 - - Nonlinear integral equation 



G , s / 1 for x > 
Step(x) = | Q for x<Q 

Eq. ( 14.121 ) is called counting equation and it relates the number of roots and 
holes appearing in a given Bethe state to the total spin of the chain. 



4.2 Nonlinear integral equation 

In the counting function one may separate the contributions of all types of 
roots and rewrite 



Z N {&) = N [<p 1/2 (tf + 9) + <p 1/2 (0 - 6)] + tp H+ /2(&) + fH-/2@) 

N H M c M w 

k=-N H k=-M c k=-M w 

N H +M R 

J2 <Pi{#-Xk) + <Pi{#) + <Pi( 2 #)- ( 4 - 13 ) 

k=-N H -M R 

In order to get a more explicit and useful expression for Eq. (14.131 ) one has 



to perform some standard manipulations [57], 58|, [59|, [60, [61, 62J] . Let notice 
that the function Z^i^d) has a very complicated analytic structure due to the 
logarithmic cuts of the functions </?!,($). To deal with a meromorphic function 
it is useful to consider the derivative Zn(i9) of the counting function. The 
calculation is based on the fact that, assuming x to be a real solution of the 
Bethe Eqs. ( 13.321 ) such that Z N (x) ^ 0, in a small complex neighbourood of 
x : fj, = x + e, it holds 

1 -ie lZN ^ + ^Z N (x + e) 

- + regular terms. 



H--X. 1 - e iZ ^( x+e ) 

For any analytic function /(x) in an appropriate domain around the real 
axis, and <£>„(i?) is, one can apply the Cauchy theorem and obtain 

* KJ / rx 2vri/i-x 7r x 2vr2 vw 1 - e^M 
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where T x is a curve enclosing x. The Bethe real roots are finite in number 
(N is a finite number), hence it is always possible to find for every root a 
closed encircling curve T x .. The sum of all these contours around points Xj 
in the real axis may be modified to a unique curve T encircling all the real 
roots and avoiding the complex ones, which is possible because in the large 
N limit the complex roots remain finite in number. Therefore 



M M 



d/i -ie iZ »MZ N (ti) 

k=-M k=-M J i 



0i (# — x + ze) 



'DC' 



2m 1 - e^^'-") 



+ °° dx . , n . . -ie iZ ^ x+ie ^Z N (x + ie) ,„„„, 
cpi (v — x — ie) .„ , , . , 4.14 

where the curve T has been supposed to be contained in the small strip around 
the real axis < e < min{7T, np, |3tn(cfc)| V k} whose parts at ±oo do not 
contribute because of the vanishing of 0i($) for very large $: this explains 
why it is possible to calculate the contour integral as a sum of two integrals 
from — oo to +oo. With the help of simple manipulations the derivative of 
the Eq. ( 14.131 ) may be rewritten as 



+00 

dx 8(x)Z N {$-x) + K{$-x)Z N {x) = 

oo 

= N [^ (tf + 6) + 0! (0 - 0)] + &(■&) + ^(20) 

N H 
+ Y, <pi(# - h k ) + <j>H + /2(#) + <j>H-/2(#) 



k=-N H 

M c M w 

} J ^i(i? - Cjfc) - Y y?i(^-Cfc 

k=-M c k=-M w 

dx 

: <pi (d — x + ie) 

-oo 
™+oo 



2m 
dx 



—ie 


iZjf(x 


- i€) Z N (x - 


ie) 


—ie 


1- 

iZff(x 


piZisr(x—ie) 

+i ^Z N (x + 


ie) 



oo 



ipxi-d-x-ie) .„ ; V, (4.15) 

2m I — e iZ N (x+ie) v > 



^ 



4.2 - - Nonlinear integral equation 



with 



K(x) = 



WW 



1 



sin 27 



27r 7r cosh 2x — cos 27 
Finally, using Fourier transform, one gets an integral equation satisfied by 
Z N (&) 



Z n {$) = N 



1 



+ 



cosh (# + 9) cosh (i? - 0) 
+ 2tt J^ G(??-/i fc ) + 27r J B(x| #+,#_) 

k=-N H 

M c M w 

-2tt g G($-c k )-27i J2 G H {d-w k ) 

k=-M c k=-M w 

N S 

-2tt J2 [G($-s k + ie) + G($-s k -ie)} 

k=-N s 

I dxGffi — x — ie)— log [l — 
ax L 



■^\M S c e iZ N (x+ie)l 



+ i [ dxGid -x + ie)— log [1 - (-i)Msc e iz N {x-ie)} u 16 x 

J dx 

where the logarithms in the last two lines are considered in their fundamental 
determination and come from the identification 



±(_l)M S c ie ±iz N (x±ie)z i x ± i 



ie) 



1 g±«Zjv(a;±*e) 

= i log [1 - (_1)A'sc e ±^(x±ie)j -p 2yri(5 ( x _ ^ 
Lit// 

with Sfc denoting the position of special roots which take into account that, 
if the function 1 — e lZ N{x+%e) cr0 sses the cut of the logarithm, the function 
log [l — e * z iv( x +* e )J has a jump by ±2ni at the exact real point corresponding 
to a special root. The following notations have been introduced as well: 



B (x \H+, H_)=F (x, H+) + F (x, #_) + J(x) + G(x) 



(4.17) 
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F(x,H) 



+00 



dk ikx . sinhf (p+1- |ff |) fc 

27r sinh ^pk cosh |fc 

for |3m(a;) 



<£w 



(4.18) 



J(x) 



+ OG 



dfc ifca . sinh \{p-l)k cosh f(p + l)fc 
27T 6 



sinh ^pk cosh |& 
for lUmfx) 



< -min(l,p) (4.19) 



G(») S / "£ 



jfc;r 



sinh |(p — 1)A; 
2 sinh |pfccosh |A; 

for |3m(x)| < 7rmin(l,p) 



(4.20) 



G//(t?) 



27rp 



coth ^SignJm(^) + coth 



zSignJm(^) 



1 



(J7r-i9)Sign3m(i?) 
P 
1 



sinhi? sinh(i9— J7rSignCIm(i9)) 



if p > 1 
if p < 1 



(4.21) 



The contribution of wide roots Gjj(il}) is different from that of close roots, 
because of a pole of G($) at $ = ±iir. To obtain the function Zjv(i?) an 
integration in -(9 is needed. If one defines 



X(tf) = 2tt / dx G{x) 



P(<d\H + ,H^)=2ir dxB(x\H + ,H_ 
Jo 
the final form of the counting function is 



(4.22) 
(4.23) 



Z N (d) = 2iVarctan 



sinhf? 



P(#\H + ,H-)+g(0\{0 k }) 



coshO 
+ i J dxG($ -x + ie) log [1 - (-l) M sc e iz N (.x-ie)^ 

- i J dxG($ -x-ie) log [1 - (-l)Msc e iZ N {x+ie)^ + u ( 424 ) 



70 



4.2 - - Nonlinear integral equation 



with uj an integration constant and 

^(^|{M)=E Cfc [X(fe)(^-^)H-X«(^H-^)] ( 4 -25) 

k 

where {■#&} is the set of positions of the various objects (holes, close roots, 
wide and special roots) characterizing a certain Bethe state, the coefficients 
Cfc are given by 

{ + 1 for holes 
— 1 otherwise 

and for any function /($) we define 

{fii($) for wide roots 

f (<& + ie) + f (■& — ie) for specials 
f(d) otherwise 

where fnifi) is defined as 

f m ! /(0)+/(0-^Sign3m(0)) if p>l 
A "" I /(#) - / (0 - ivrpSignJm(^)) if p < 1 

Note that the previous definition ( 14.211 ) of Gn{d) is in agreement with this 



more general one. From the analysis given in the periodic case [57], ]6_ll, 163] , it 
follows that the number Nh — 2Ns remains constant along the volume size, 
therefore, one can refer to this number as the effective number of holes N e ff. 
Notice that for real x and e — ► 



Q N (x) = lim2Jmlog [l - e iZN(x+ie) ] = 

1 / 1 — e iZ N (x+ie) \ 

= lim - .-, ._,_., = Z N (x) + vrmod 2n . (4.26) 

It is possible to compute exactly the integration constant u as well. In fact, 
one has to take into account the asymptotic values of functions x($) 

Xoo = x (+oo) = ILl^l (4.27) 

2 7T — 7 
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and the asymptotics of the boundary term P {d \H+, H- 



P^P( +0 c\H + ,H_) = *- i *-^ H + + H - i \ (4.28) 

I 71 — 7 

This allows to calculate the aymptotic value of Zn($) and, making use of the 
counting equation ( 14.121 ), it follows 



w = (ir- 2 7 )S + 4tt - 7 (H+ + H^ + 4) + vrSign(j9 - 1)M, 



u 



- 2 Xoc S - Poo -A Xoo + ^ 1 (H + + fl_) + 2 -^lL(S + 2) = 

7T 7T 

= 7rSign(p-l)MH, . (4.29) 



The last term is a multiple of 2ir if M^ is even, which is the case unless 
there are an odd number of self-conjugated roots. Zjq($) does not change if 
one adds or subtracts 2tt because of the quantization rule. Without loss of 
generality one can take 



to = n (M sc mod 2) (4.30) 

therefore, it is convenient to include this constant into the logarithmic term of 
the convolution, by modifying it to log [l — (— l) (S e jZjv *- x '] , with 5 = Msc mod 2. 
The quantization rule of any state becomes 



Z N (&j) = 2ttI 3 , /,-6Z+-. (4.31) 

Similar considerations can be done to express the energy in terms of the 
counting function and the positions of the various roots. The expression 
which one derives depends on the number of lattice sites N and the lattice 
spacing a, it has the form 
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2aE^ = Y J C k { 



1 



1 



cosh(e-tf fc ) cosh(e + $ fc ; 



dx 

2^ 



sinh (G — x) sinh (— G 



x) 



cosh 2 (G — x 

_ w yg. (e ,,)( 



cosh 2 
1 



e-x) 



Qn{x) 



dx 



cosh (Q — x) cosh (Q + x) 

2tt5(x)] 



where 



—$(e,x)[P(x\H + ,H_ 
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$ (0, x) = <p 1/2 (6 - x) + 01/2 (6 + x) 



(4.32) 



4.2.1 Continuum limit 

At this point one can calculate the continuum limit of the energy, i.e. an 
expression valid for iV — > +oo, a — > in such a way that the finite volume 
L = Na remains fixed. One can try to do that by computing first the 
continuum limit of the counting function 



Z{fl) 



lim Z N U)). 



One has to set 



O ~ log 



2N 
ML 



(4.33) 



(4.34) 



with M. a mass scale. The source term g{'d\{ q &k}) and the convolution 
integral do not present any dependence on N, which is instead crucial in the 
driving term. As a result, one can define the Nonlinear Integral Equation 
(NLIE) 



Z(&) = 2ML sinh ■& + g (•& \{& k } ) + P (•& \H+, H_ ) 

- 20m f dxG(tf -x-ie) log [l - (_l) V z(x+fc) ] 



(4.35) 
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For the vacuum state containing real roots only, this equation coincides with 
the one found some years ago in [58|. The quantization relation 

Z{0 i ) = 2*I j , I j eZ+ S - (4.36) 

holds as well and fixes the positions of holes and complex roots determining 
the source term g ($ | {$&}) characterizing the Bethe excited states. The term 
P (x \H + , H_ ) takes into account the effects of boundaries and, in particular, 
it depends on the parameters H± which can be related to the DSG boundary 
parameters <fi±. We will be considering this point in the next Chapter. 

In order to compute the continuum limit of the expression ( 14.321 ) for the 
energy, let us just say that the following asymptotic behaviours hold for 
N -> +00 

— arctansinh (0 — •&) ~ — Me + ^ 

LJ LJ 



N f dx Qn(x) M f dx _ 



LJ 2tt cosh (±6 - x) 2 J 2n 6 ^ (x) 

The last formula is the contribution of the integral term in the second line of 
Eq. ( 14.321 ) with Q(x) = lmiTv^oo Qn{%)- The contribution of the third line 
of Eq. ( 14.321 ) can be evaluated by considering that it takes the form 

<N) _ 2N 2 f dx d ( y?i/2(» 



E ^-~ L J tt d6 Vcosh(26-x)J- (4 ' 37) 

It is, evidently, divergent for N — > +00. The contribution of the boundaries 
to the energy is given by the last line of ( 14.321 ) and can be expressed, making 
use of Fourier representation, as 



(jv) 2iV f dk coskQ 



J B 



L ./ 27r cosh |& 



sinh \{j)+l-H + )k 
sinh I (j> + l)fc 



sinh I (p + 1 — H_ ) k sinh | (p — l)k 
sinhf(p + l)A; sinhf(p + l)A; 



(4.38) 



This term, after the shift k — > ik, splits into two terms. The first one diverges 
with N, let us call it E^^. The second one takes the form 
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M (sin^ip + l-H^ 
Eb = 



2 V sinf(p+l 



sin f (p + 1 - #_) ^vr 



cot T (p+l)-l (4.39) 



sinf(p + l) 4 

At this point it is allowed to renormalize the energy subtracting all the di- 
vergent ' 
result is 



vergent terms, i.e. shifting the energy E —* E — E B nlk — i^uik- The fi na l 



E = E B + M J2 Ck cosh (fc) # k - M J ^ sinh(x) Q(x) (4.40) 



Eqs. ( 4.351 ) and ( 14.401 ) are the most important results of this Chapter. In- 
deed, they allows to express in a continuum framework the energy, i.e. the 
eigenvalues of the Hamiltonian ( 13.301 ) in the formalism of Chapter 3, as a 
function of Bethe roots and boundary parameters. The energy ( 14.401 ) is the 
perfect candidate to define the scaling functions of the model and to relate 
different states in the IR and UV limits introduced in the next Chapter. 
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The continuum theory 



Once the Z($) and the energy are known, in the IR limit A4L — > +00 one 
can make contact between the NLIE excited states and the scattering theory 
of the underlying field theory. On the other hand, in the UV limit ML — > 0, 
one finds contact with the CFT at the UV regime of the DSG model. 



5.1 IR limit 

The IR limit is expected to reproduce a field theory in a large volume L. As 
a scale parameter it is convenient to choose the dimensionless size parameter 
I = M.L with Ai playing the role of a mass. Therefore, the very large 
I regime can be interpreted both as infinite large size, that reproduces the 
conditions of a field theory in an infinite Minkowski space-time, or large mass 
scale, compared to energies in play, that corresponds to an infra-red point 
in the RG flow. There are many possibilities to give an interpretation of the 
physical content of the model described by the NLIE. The main question 
is to find a correspondence between the physics of the IR model and the 
different states of the NLIE which one can deal with. Let us notice that the 
infinite volume limit is not a thermodynamic limit, because the number of 
the various objects in the NLIE, apart the roots on the real axis, remains 
fixed. So, it is expected that the energy E should approach a finite limit 
depending only on a finite number of roots, reproducing thus a free massive 
spectrum. In order to prove this conjecture, it is sufficient to notice that for 
I large enough the integrals in Eqs. ( J4.351 ) an d ( 14.401 ) become small of order 
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O (e l ). They can be dropped. Hence 



Z(#)=2lsmh'& + g(&\{& k }) + P(#\H + ,H_) for I -> +00 (5.1) 

which can be seen as a constraint on the asymptotic states. Similarly, in the 
energy formula the integral term drops and 

E = E B + M^c k cos\ k) $ h (5.2) 



5.1.1 Holes- Solitons correspondence 

The important thing is to understand what type of constraints Eqs. ( 15.11 ) and 
( 15.21 ) fix and what they mean from a physical point of view. Consider the state 
without complex roots [Mc — Mw — 0). Thus one can take into account two 
different situations. The simplest one is that corresponding to the presence 
of real roots only. It is evident that in this case the energy reduces to E = Eb 
because neither the NLIE nor the energy depend explicitly on these roots, 
i.e. one has a state with a fixed energy value in which from the constraint 
( j5.ll ) it does not come any informations. One could associate this NLIE 
state to the vacuum of the theory. This assumption is also supported by the 
fact that for / — > +oo the number of real roots goes to infinity forming a 
vacuum Dirac sea and that, from the Eq. ( 14.121 ), the total spin of the system 
is S — 0. This is exaclty true only if \h±\ < h c where the quantity h c has 
been defined in Chapter 3, otherwise the behaviour of energy presents some 
special effects and, in order to have S = 0, one should consider the presence 
of some different type of roots. We will be treating this situation later, so 
far let us consider only the range \h±\ < h c . 

The other possibility is to have a NLIE state defined by Nh holes in 
the positive real axis, say at positions $&, k = 1,2, ... Nh- For sim- 
plicity take just two holes at positions $i and i?_j. Due to the fact that 
{$-k} = { — $k} V k = 1,2, ... and to the antisymmetric property of the 
NLIE, for any hole of type $_£,. the same equation holds and one can con- 
sider only objects with positive real part: each of them implicitly describes 
also the contribution of those with negative real part. From now on we 
will consider objects with positive real part and use the notation Na with 
A = H,C,W, ... to refer to their numbers. The NLIE valued at $i becomes 
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Z (0i) = 21 sinh i?i + x (2^i) + P (tf a |#+, #_ ) 

= 2Zsinhtf 1 +.F(tf 1 ,i/ + )+; zr (tfi,i7_) = 2tt/i (5.3) 



where 



1) 



hoo 1 j ™„V, 37T ; 



F(x,tf)4 / d^sinhf&sinhftp-l^ 



(5.4) 
Eq. ( 15.31 ) may be interpreted as a quantization rule for the momentum of a 
particle of energy Ai coshi?i. In fact, taking the exponent 

e 2iZsinh#i e i[:F(tf 1 ,H + )+.F(i?i,ff_)] = ]_ ^ gN 

after some mathematical manipulations [64], one can prove that 

F($,H) = -ilogK + ('& : H) (5.6) 

where TZ + ($,H), which has the same structure of the function V+(u) in- 
troduced in Eq. ( I3.14I) , is exactly the soliton boundary reflection matrix as 



found in [24 1: therefore 



e 2iisinho in + ^ H+) n+ ^ H j = j ^ ^ 

which corresponds to the quantization rule of a particle moving in a volume L 
with momentum M. sinh $i towards two boundaries and interacting off them 
with reflection amplitudes 1Z + ($1, H + ) and 1Z + ($1, HJ). 

The parameters of reflection matrices are related to the Dirichlet bound- 
aries (f)± in the boundary sine-Gordon Action and the boundary magnetic 
fields of the XXZ based construction, then, taking into account that <f>± are 
always positive parameters while the parameters H± and h± take values in 
positive and negative ranges, one gets 

P (i-U ± ) if h ± >o 

H± = { > P I 5.8 

P (l + |0±) if h±<0 

In particular, from our choice of boundary magnetic fields (/i_ > and 
h+ < 0), it follows #_ = p(l- |0_) and H + = p ( 1 + |0 A It is thus 
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tempting to say that the one hole NLIE state corresponds to a soliton of 
DSG. This interpretation is also supported by the fact that the energy of 
this state is 

E = E B + Mcosh$ 1 (5.9) 

that looks like a particle excitation over the vacuum (remember the parametriza- 
tion of particle energy in terms of rapidity given in Chapter 2). Moreover, 
from the counting equation it follows that the state has spin S — |, i.e. 
an excitation over the vacuum spin S = 0. We conclude that the physical 
meaning of a hole state is a soliton with energy M. cosh $i scattering off the 
boundaries via the Ghoshal-Zamolodchikov reflection amplitudes |65f . 

Consider, now, a more complicated state, say the one with a generic 
number Nh of holes at positions $k, k = 1,2, ... Nh- Due to ( 13.91 ) and to 
the definitions ( 14.201 ) and ( 14.221 ), the function x(-#) can be written as 

x (tf) = -iloga(tf) = -z log £++(#) (5.10) 

where £++($) is the sine-Gordon soliton-soliton scattering amplitude. There- 
fore, for Nh holes one has that the spin of the state is S = —^-, the energy 
is 

N H 

E = E B + M^2cosh$ k (5.11) 

k 

and for any hole dk the following equation holds 



N H 
N H 

= 21 sinh § k + J2 [X 0?fc - 0j) + X (^ + $i)\ +F(#k, H+) + T{d k , #_ 



(5.12) 



After an exponentiation, Eq. ( 15.121 ) becomes 



N H 



]JS ($ k \# j )K + (# k ,H + )1l + (# k ,H_) = 1 (5.13) 
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where, for simplicity, we have defined 

S (0 \dj ) = S++ (■& - ^) S+X (0 + fy) (5.14) 

Eq. ( 15.131 ) is the quantization rule of a soliton interacting on two boundaries 
via 1Z + ($£, H + ) and TZ + (i?^, if_) and with iV# — 1 solitons via the soliton 
S matrix S++ ($k ± i?j). Notice that the signs ± refer to the fact that the 
soliton #fe can scatter with the other solitons either before or after they have 
respectively scattered off one of the two boundaries. 

5.1.2 The antisoliton 

According to the symmetry C introduced in Chapter 3 we can deal with 
states of negative spin just reversing both magnetic fields h±. In terms of 
DSG model this symmetry requires some care. It corresponds to reversing 
the topological charge Q — > —Q, but, in order to achieve the correct result 
it is not enough to send <fi± —>■ —<f)±. In fact, an exact definition of the 
DSG theory describing a state of topological charge — Q in the same ranges 
of boundary fields, comes from the successive application of the symmetries 
— > —(f), <p — ► + 7f and <f>± — > — (f>±, (f>± — > (fi± + ^ simultaneously because 
one has fixed the boundary fields to be always positive. This leads to the 
same description in terms of NLIE states for physical states with charge — Q: 
a state with Nh holes describes iV^ particles of mass M. scattering with 
S-matrix given by the 5 , ++(^) = SZZ("d) sine-Gordon amplitude. Therefore, 
one espects to get a theory describing corresponding conjugated states of 
those treated previously, for instance antisoliton states instead of soliton 
states. From Eqs. ( 15.81 ) it follows that, due to the transformations of fields, 
in all formulae one has to substitute 



H± + l-*-H±-l , V H± : -p- 1 < H± <p+l . (5.15) 

As a consequence, the soliton reflection factor goes into the antisoliton reflec- 



tion factor 1Z~ (i!},H±) found in [24] . Under these transformations the pure 
holes state can be interpreted as a state of Nh antisolitons of sine-Gordon 
with Dirichlet boundaries d>+ 



D 2il sinhi?^ 



l[S($ k \^)n~ (& k ,H + )K- (#*,#_) = 1 (5.16) 
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where S (-d k \&j ) has been defined in ( 15.141 ). 

5.1.3 Boundary bound states 

From the analysis given in [66] of the Bethe Ansatz equations there is evi- 
dence that boundary bound states should be described by purely imaginary 
objects. Therefore, in the continuum NLIE we search for solutions charac- 
terized by the presence of one or more roots on the imaginary axis. Such 
roots give contributions X(k) ("$ — ^ a k) to the NLIE source term g ($ \{i9 k }) if 
their position is at {} k = ia k . Moreover, in the large I limit from the oddity 
of the Z($) it follows that any such a root implies 



'ftzZ (ia k ) = 2nl ak = 
JxnZ (ia k ) = 21 sin a k + 3mx (2ia k ) + 3mP [ia k \H + , f/_ ) = 

The exploding term 21 sin a k for / — > +oo has to be compensated by some sin- 
gularity coming from the other two terms. This may happen when one consid- 
ers poles in the boundary reflection amplitude. In the function P (i? \H + , H_ ) 
the boundary content is encoded in functions F($,H±), hence, only poles 
coming from these terms and depending on boundary parameters H± can be 
taken into account, any other poles has to be interpreted in the bulk theory. 
In the region < H± < p + 1 there are no poles in functions F (ft, H±). All 
possible poles in the integral terms of the NLIE may be interpreted in the 
bulk theory. Hence, taking into account the range of definition of boundary 
magnetic fields h± which we have chosen, the poles which one is interested in 
have to be found in the region — 1 < H^ < for the boundary parameter H_ 
and — p — 1 < H + < — 1 for H + , this means that /i„ > h c and h + < —h c and 
clarifies why in the previous section, where \h±\ < h c , we have not treated 
boundary bound states. To find such poles one has to consider the large |A;| 
limit in the corresponding Fourier transform F (k, H). It follows 

F(k,H)~ e -% mkl for |£;|^+oo (5.18) 

therefore, in F ($, H) there appear singularities at 

9te(tf) = , 3m(tf) = ±- (H + 2pl) with lei. (5.19) 
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These singularities correspond exactly to the first set of poles ( 13.181 ) parametrized 
by Dorey and Mattsson in [48] if H± take values in the range ( 15.81 ): — 1 < 
H_ < and — p — 1 < H + < — 1. The sets of these poles can be denoted as 

A ± = jz«W = _j| [H±±p + p(2n + 1)] : 

. >0, n = 0, 1, . . . ^n*} (5.20) 



7T 

2> 



a?" 



> 



ai^ 



where n are integers defined such that < \H± ±p + p (2n + 1)| < 1 



and the roots z<% = i\ [H± ±p + p(2n + 1)] may be discarded because of 
the properties of the Bethe equations. So one can consider the boundary 
bound states as states in which combinations of purely imaginary roots with 
imaginary part in the set ( 15.201 ) appear. 

Let us deal with the simplest case with only one excited boundary, say 
the one at the left. In terms of NLIE this means the presence of a pure 
imaginary root ia ~ = —i^H^ which fixes the expression of the IR energy 
to be 

E = E B -M cos a { ' ] (5.21) 

where, due to range of definition of the boundary parameter H_, the term 
M. cosOq is always positive, i.e. the energy becomes smaller than the vac- 
uum energy Eb ■ In general, the whole set of pure imaginary roots contributes 
a negative amount in the energy and these roots cannot be considered as ex- 
citations. This means that the vacuum which we have been working with 
is instable for h- > h c and h + < —h c . The vacuum decreases as the value 
of boundary parameters take values grater than the critical one. This could 
appear as a non-sense, but remember that the vacuum energy has been de- 
fined up to a bulk term, therefore, in principle, one could redefine that by 
shifting by a finite amount. In this sense any root on the imaginary axis 
corresponding to a pole of F (■#, H±) defines lower levels of vacuum energy. 
To cure this situation one has to define a new ground state by adding up to 
the Dirac sea made of real roots all the pure imaginary roots in the set A ± 
defined in ( 15.201 ). The vacuum energy becomes 



E B bs = e b -mJ2 cos 4 +) - m J2 cos a i~ ] > 4 ±} e A± ( 5 - 22 ) 



fe=0 fc=0 
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For reasons of consistency, one expects the total spin associated to the vac- 
uum state to be S = 0. This is possible only by redefining the formula ( 14.31 ) 
as S = y ~ M = y ~ Mr ~ Mi, where Mr is the number of real roots, as 
usual, and Mi the number of pure imaginary roots in the set A . This new 
definition of the total spin might have some consequences on the counting 
equation. Indeed, it has been derived taking into account the relation among 
the spin, the number of lattice sites and the number of total roots. Notice 
that from the first of Eqs. ( 15.171 ) it follows that the quantization numbers 
of pure imaginary roots are always zero. This means that the definition of 
-fmax, the greatest integer smaller than Z N {+oo)/2n, and 7 min , the smallest 
integer greater than Zn(—oc)/27t, needs some care. The quantization num- 
bers of pure imaginary roots, due to the fact that they are zero, do not give 
any amount neither to the sum of all quantum numbers nor to the greatest 
or the smallest among them. Therefore, i max and i m i n need redefining as 
I max = 2S + M -Mi, I min = -2S -M + Mj. Following the same deduction, 
the number Ni of pure imaginary holes does not depend on the integers i max 
and imin, hence i max — imin = 2 (Mr + Nh — 2Ns) as in ( 14.111 ). As a result, 
the counting equation remains the same and the numbers of pure imagi- 
nary roots and holes are not involved in it. In particular, one can choose 
their number independently from the total spin and the numbers of all other 
objects. 

In this scheme boundary excitations are obtained by successive removing 
of such imaginary roots, i.e. by creating holes ■$& with d\t (■#&) = and 
i(3m('&k)) G ^ ± - The creation of a hole at position, say, ia = —i\H- 
describes the excited boundary with energy coming from the Eq. ( 15.21 ) 



E = E B bs + M cos a { ~ ] (5.23) 

This is the first boundary bound state at the boundary on the left. Notice 
that it is exactly an excited state over the new vacuum E B hs . In particular, 
the state described by a number of imaginary holes equivalent to the num- 
ber of imaginary roots continues to be an excited state although its energy 
becomes quantitatively equivalent to E B , in fact, for \h±\ > h c the system 
is completely different from the case \h±\ < h c and the quantity Er has not 
longer the same meaning. 

One also expects that a boundary bound state implies some consequence 
on the scattering theory, in the sense that particles scattering off the excited 
boundary have to be affected by a different behaviour from that of free par- 
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tides in not excited boundary system. Let us, thus, consider the simple state 
consisting of an excited boundary (an imaginary hole at position — i\H- and 
a real hole at position $i). Substituting in the NLIE, considering its antisym- 
metric properties and the quantization equation, it follows that the position 
$i of the hole fulfils a quantization rule of a particle of energy A4 cosh^! in 
a finite domain such that 

e 2,lsmMl ? ( + _ 0) fa, HJ) TZ + {A, H+) = 1 (5.24) 

where, making use of the function ( 15.141 ), the notation 

V+ ± )0) (■&, H) = S (■& \ia ( ^ ) TZ + (0, H) (5.25) 

has been introduced with a G A^. It corresponds to the first boundary 



bound state reflection amplitude given in [48|. From a physical point of view 
one can look at this result as a soliton in the boundary system scattering on 
the left excited boundary by the amplitude Vt_ Q \ ($,H_) and on the right 
one by the unexcited amplitude 1Z + ($, H + ). The energy of this state turns 
out to be 

E = E% bs + M cos a { ~ ] + M cosh A . (5.26) 

From now on we will refer to the amplitude Vt ± -, (i?, H) as the solitonic first 
excited boundary reflection factor. 

Boundary bound states may occour at the right boundary simultaneously 
as well. In this case, evidently, a first pole in F ($, H + ) may be found at po- 
sition ia = — z| (H + + 2p). The two excited boundaries state is described 
by creating two pure imaginary holes at the exact position of the two imag- 
inary roots i«Q , ia and 

e 2UM V+_ fi) (A, H.) P ( + +i0) (A, H+) = 1 (5.27) 

with energy 

E = Ef* + M cos a { ~ ] + M cos a[, +) + M cosh A . (5.28) 

In this case the reflection amplitude at the right boundary is excited as well. 
Obviously, when \h±\ > h c any real hole can be regarded as in the previous 
section such that Njj of them describe Nh solitons reflected back of the 
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boundaries and scattering each others by the sine- Gordon S matrix. If the 
boundaries are unexcited they are reflected by the usual soliton boundary 
amplitude TZ^ (•&, H±) through the Eq. ( 15.71 ) and a reference ground state 
energy E^ s , otherwise by V?± s ($, H±) from the first excited boundary. 
Considering such a state, it is simple to find that, after an exponentiation, 
the NLIE becomes for any real hole 



N H 

e 2ii s ^ k "Q s ^ |^. } v +_ q) ^ H) v + + Q) ^ H+) = 1 , (5.29) 

According to the values of H± and p, it can happen to have more than one 
pure imaginary hole on the left and on the right boundaries. In this case, 
considering a number n of pure imaginary holes in the set A - and n' in the 
set A + , indeed these numbers are in general different because of the strong 
dependence of positions of pure imaginary objects on parameters H±, any 
soliton with energy M. cosh -d k should scatter by the two excited boundaries 
via the equation 



N H 

e 2il ^^Y[S^ k \^)V+_ >n) ^ k ,H_)V^ nl) ^ k ,H + ) = 1 (5.30) 

where the excited boundary reflection factors Vt n x ($&, H_) and Vt + n ,\ (i?fe, H + ) 
on the left and on the right have the form 

Q 

V+ ±>q) {4,H±) = U S {^ \i°i ±) )K + &H ± ) (5.31) 

i=0 

and the terms S ($& |i?j ) take into account the presence of other solitons (real 
holes) . 

In order to study the antisolitonic excited boundary reflection factors one 
has to consider the symmetry which changes states of topological charge Q 
into states of charge —Q. It has been stated that this corresponds to send 
H± + 1 — ► —H± — 1 with the consequences that 1Z + (t?, H) — > TZ~ (i?, H) and 
one deals with antisolitonic states. Observe that the antisolitonic reflection 
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factor 1Z (■#, H) has simple poles at position ia_^ as the solitonic reflection 
amplitude. In addition it presents simple poles at positions 



7T 



SKe(tf) = , 3m(0) = T~ (2 + H + 2pj) with j e Z (5.32) 

which represent a second set characterizing new boundary bound states 



*fl 



(±) 



7T 



i- [2 + Jf ± ± p - p(2m - 1)] : 



TV 

2 > 



fl 



(±) 



> 



# 



(±; 



... > 0, m — 1, 2, . . . . < m 



± 



(5.33) 



where m^ 1 are integers defined such that < |2 + H± ±p — p (2m — 1)| < 1. 
According to [48] these poles represent the second set of boundary bound 
states ( 13.191 ) describing antisolitonic boundary bound states. Therefore, one 
can repeat the same derivation given above for solitonic states and find similar 
equations. The most general case with Nh antisolitons and boundary bound 
states at left and right boundaries characterized by m and m' pure imaginary 
holes in sets B can be expressed through the following expressions for the 
quantized scattering of a particle of energy A4 cosh dk and for the total energy 



N H 



J2il sinh$ ; 



J] S {<& k \d 3 ) P(l m) (i) k , HJ) V^ m>) (0 fc , H+) = 1 (5.34) 



&k 



N, 



E = El hs + MJ2 cos A- +) + M J2 cos A H + M J2 cosh #j ( 5 - 35 ) 

j=l i=l j=l 

where 



r 

P ( - ±r) (0, i7 ± ) = J] S (f \i/3F> ) K- (0, H ± ) 

j=0 



(5.36) 



with r = m for the left reflection amplitude and r — m' for the right reflection 
amplitude. The reference ground state is redefined as 



El hs = e b - M JT cos A (+) - M J2 cos A 



(-) 



(5.37) 



i=\ 
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5.2 UV limit 

Now we want to give a computation of the I — > limit in order to describe 
the system in the UV limit. Starting points are the NLIE and the expression 
for the energy ( 14.401 ). 

It is known from boundary perturbation theory that the scaling behaviour 
of the energy for small volume size is 



E(L)~l(*-±) (5.3* 



as sketched in Chapter 1 and Chapter 2. In this limit the roots and holes 
may rescale to infinity or stay in a finite region. Their behaviour is 

tf fc = 0± ± log - or d k = #1 (5.39) 

where thus we classify these objects into three types denoted by the indices 

"±, 0" 

K(0} = {^±logy} , {<(/)} = K}. (5.40) 

Accordingly the NLIE splits into left and right kink equations. Define, thus, 
the kink NLIE as 



Z±(t?) =limZ(^±logy 
and introduce the scaling function 



Z (ti) =\imZ(ti\l) (5.41) 



Ei(l) = ^boundary ~ 7^7^(0 wlth }™ C *W =C ~ 24A ' ( 5 " 42 ) 

where -^boundary refers to any contribution coming from boundaries and not 
depending on the space size L, that is in the limit I -^ the boundary bulk 
energy terms can be ignored because they are finite quantities. 

Due to the antisymmetry of the NLIE it is enough to consider just the 
right kink equation and the objects with positive real parts. All the others 
have to be described by the same equations and one gets for their contribu- 
tions the same results. Observe that the source terms corresponding to the 
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various roots and holes may have different behaviour according to their set 
of reference ( 15.401 ). For source in {i?j!~(0} 

limx (# + log y - ^ + (/)) = X (0 - tfjJO (5.43) 

and for sources in {$*!(/)} 

lim X ^ + logy-^(0J =xf^ + logy-^J =Xoo (5.44) 

where Xoo has been defined in ( 14.271 ). Therefore, given e.g. Ng holes of 
type "+", there will be in the NLIE Nh — N^ contributions of term Xoo- 
One can introduce for the numbers of objects in the sets ( 15.401 ) the notation 
Np , Mp , . . . and expect that counting equations similar to that given in 
( 14.121 ) hold for the spin S + and S° according to the condition that Ng = 
2N+ + 2N° H , M c = 2M+ + 2M° , . . . and S = 2S + + 25°. Here the factor 2 
take into account that the number of objects with negative real part is equal 
to that of objects with positive real part and that each of them carries the 
same spin. Hence, one considers only sources with positive real part and take 
care about the factor 2 to have the total contribution of any source. Notice 
that there might be Nj contributions Xoo as well, due to pure imaginary 
holes. It has been stated that they do not count in the counting equation 
and cannot be related to the spin, but they have to be taken into account. 
After few manipulations one has, for I —* 0, 

Z + (d) = e +d + g + ($\{$+})+P(tf\H + ,H_)+ f dxG($ -x)Q + (x) (5.45) 
where 



g + (#\{&+}) =KmgU + hg± {0 fc (Z)} 

f N \ M+ 

= 2 Xoo (s + ^-S + )+J2 CkX(k) (0 - 0+) + vrL+ (5.46) 
^ ' k 
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where L^ = Sign(vr - 27) (M w - M&) and M+ = N+ + M+ + M+ + iV+. 
Similarly for sources that remain in a finite region 

Z o {0) = g (# |K}) +P($ \H+,H.)+fdxG{-& - x)Q {x) (5.47) 
where 



g ($\{#° k })=\img(#\{Ml)}) = 

= 2 Xoo (S - S°) + J2 CfcX(fc) (^ - **) (5-48) 

k 

The functions Q +t o(x) are related to Z+ t o(x) as Z(x) to Q(x), M° = N^ + 
M% + M^r + N% + JV>. In order to calculate the UV of the energy ( 1440T ). 
some asymptotics for i? — > +00 are required: 

Z + (+oo) = +00 , Q+(+oo) = (5.49) 

g + (-00 J {tf+}) =2 Xo o(s + ^-- 2sA + 2ttK+ (5.50) 

Z+(-oo) = <?+ (-00 I {^} ) + Poo + — Q+(-oo) (5.51) 

7T 

with i^ e Z and P M defined in formula ( |4.28j ). By definition (see Eq. 
(S255), for any K E Z, it holds 

Z + (±oo) = g + (±oo) + 7r + 27rir (5.52) 

therefore, making use of the asymptotics ( 15.491 ) and ( 15.501 ) and comparing 
Eqs. (ME) and (5321 ). 

q^-oo) = * ^ < S + J - 2S+ > + P " - T - 2 ^ + ^ (5.53) 

and finally, observing that for / — > the terms in the energy behave as 
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M. f dx 1 f dx 

T y-sinhxg(x)~ z y-e-g + (x) 



M. cosh dk ~ 
one can rearrange the formula ( 15.381 ) as 



L e +K 
L ' 



c 1 

A = 1 

24 2vr 



'M+ 
. k 



^ ~ j t e+XQ + (x) + ^£ COS ^ (5 ' 54) 



where the last term takes into account the whole set of pure imaginary holes 
i'dk in the sets A and B ± describing in the IR limit the boundary bound 
states. For / — > this term goes to zero. At this point, making use of the 
notations introduced above, defining the sums over all quantum numbers for 
any type of sources as I\ = V • l\- with A = H,C,W, ... and following 
the scheme given in |57j, one gets the following expression 



V 



p+1 
-2 Is 



where 



2p 



1+ 



V 



2p 



S 



Nj_ 
2 



+ N 



(5.55) 



N = J+ -/+-/+- 2J+ + L+ S + + 2k-S + -2 (S + ) 2 e Z + 

and k = K — K^. This is an exact computation of the scaling limit of the 
scaling function ( 15.421 ). 



5.2.1 UV limit as a BCFT 

The UV limit allows one to make contact with a c = 1 boundary CFT. The 
scaling function as defined in ( 15.421 ) takes the form 



1 - 24A, 



(5.56) 
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and, in general, for a theory with non trivial boundaries Aj 7^ for the 
conformal vacuum too. In our case the main aim is to relate NLIE states 
of the IR limit to conformal states. In the sections devoted to the large 
volume limit, we have established that, for \h±\ < h c , the vacuum state is 
characterized by the NLIE state without any source and with spin S = 0. In 
the language of UV limit this corresponds to S = S + = S° = 0, therefore 

It is the effective central charge in the UV fixed point, in particular q < 1. 
Eq. ( 15.571 ) describes, thus, the scaling of the vacuum of the theory and 
reproduces exactly the formula given some years ago in [66] . These few 
considerations confirm that the NLIE allows to control the complete flow of 
the vacuum scaling functions from the IR to the UV limits. To be more 
precise, if the vacuum of the scattering theory is fixed to be Eb as in ( 14.391 ), 
which is a quantity strongly dependent on the boundary parameters H±, 
the correspondent conformal state has an energy Ei = — ^f^cf 10 according 
to Eq. ( 15.421 ). It is evident that the boundaries contribute some energy to 
the Casimir effect due to the finite size. The complete flow of the vacuum 
energy E (l) from zero to infinity can be seen in Fig. ( ]5.11 ). In Figs. ( |5.2| ) and 
( J5.31 ) the behaviour of the corresponding scaling function c (/) is displaied. In 
order to obtain such pictorial representations, the NLIE has been integrated 
numerically at several values of I. 

In a more general situation, one would be able to treat the excited states 
as well. For semplicity let us introduce the notations 

R =\IW = ir • ™ = 25 + iV ' eZ > 

where S denotes the total spin carried by sources at positive real position, 
i.e. the total spin of the system is exactly twice. The conformal dimensions 
of any conformal state can be parametrized in terms of relevant quantities of 
NLIE states as 



2 



S° + k 1 



-1 2 



H — raR 



+ N (5.51 



% R 2 

where <f>± are the boundary values of DSG related to the parameters H± as in 
Eq. ( 15.81 ). Eq. ( 15.581 ) (see [65] ) reproduces exactly the conformal dimensions 
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Figure 5.1: The behaviour of the vacuum energy as a function of the scale 
parameter I for the case: p = | , H_ = 0.3642, H + = 0.4238. 
According to Eq. ( 15.42} ), Eq can be discarded. Therefore 
E{1) ^ forl^ +oo. 



in the partition function ( 11.731 ) for all states with S° + k = 0. We will see 
that any physical state satisfies this condition. Therefore, without loss of 
generality, we can rewrite for physical states 



2 



n 2 



- H — mR 

ix 2 



N 



(5.59) 



This result can be exactly compared with the c = 1 CFT of a boson with 
Dirichlet conformal boundary conditions which is compactified on a circle of 
radius R [671, [68 . 



5.2.2 Classification of BCFT states 

The Hilbert space 7i of the c = 1 CFT (free boson) with Dirichlet boundary 
conditions is composed of Heisenberg algebra representations 7i m of c = 1 
CFT whose primary states \m) are created from the vacuum by radial ordered 
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Figure 5.2: The scaling function c(l) for the vacuum state with: p 
H^ = 0.3642 and H + = 0.4238. 




Figure 5.3: Particular of Fig. ( 15.21 ) . According to analytic formulae, one 
expects an UV central charge c = 0.999618 . . ., here it is 
reproduced with a precision 10~ 6 . 
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vertex operators e %K ^ of U(l) charge and corresponding conformal dimensions 
1 /</>+ -(f)- 1 



n 2 



mR\ , A m = 2K 2 . (5.60) 

The whole set of states of a representation TC m is created by applying repeat- 
edly the creation operators ( 11.591 ) in such a way that 

H m = {a- kl . . . a- kq \m) \k\, . . . , k q e Z+ } (5.61) 

Notice that only creation operators in the holomorphic sector need applying. 
The Hilbert space is, thus, a linear combination Ti. = ® m A/A m H m of sets 
( 15.611 ). For any state \i) e TC m the energy is given by 



E i =*\A m + N i -±) , N i = J2k J eZ + (5.62) 



L 



It is evident that the value in the expression ( J5.59I ) is a good candidate to 
classify several states in relation to NLIE states, one has just to identify 
A = A m + Ni. For instance the conformal vacuum takes the form 



A vac = ^[4>+- <M 2 = A • (5.63) 

Solitonic primary states 

The first interesting conformal primary states which one can consider are 
those correponding to solitonic IR states. It has been shown that sine-Gordon 
solitons can be represented in the scattering theory limit as pure NLIE holes 
states. Take e.g. the state of one hole at generic position i?i in the positive 
real axis and take the minimal quantum number I x — 1. In the IR theory 
it describes a soliton of energy M. cosh $ 1 scattering off the left and right 
boundaries. Suppose that when / — > the hole becomes an object in the 
set {#+(/)}, i.e. its position rescales to infinity. The spin is S — | and, 
in particular, the spin has to be S = S + = \ and the quantum number 
Ii — I\ — 1. Hence, the associated conformal state in the UV limit is 



2 



l R 



n 2 



Ai (5.64) 
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Figure 5.4: One soliton state energy as a function of I. p = |, i?_ = 
0.3642, H + = 0.4238. For I — ► +oo t/ie energy goes to a posi- 
tive finite value corresponding to the soliton rest energy. For 
simplicity it has been set A4 = 1. TTie filled line represents 
the case with I\ = 1, the dotted line I\ = 2 and the dashed 
line I\ = 3. 



One is thus dealing with the conformal primary state of the Heisenberg rep- 
resentation Tii. Observe that in this situation one is forced to choose I± = 1. 
Any other choice should not be compatible with a primary state. In fact, 
any quantum number greater than 1 should imply N > 0, i.e. a descendant 
state with energy grater than E\. Therefore, the quantum number is crucial 
to have exactly a primary state. The behaviour of the one soliton energy is 
represented in Fig. ( 15.41 ) for cases with I\ — 1, I\ — 2 and I\ = 3. 

It is possible to take more complicated states. For instance the state with 
a generic number Nh = iVj, that is a state with Nh holes in the positive 
real axis lying in the set {$t} for / —> 0. According to our notations 



S = S~ 



Nh 
2 



m = N H and N = I H - -N H (N H + 1) 
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therefore 



N, 



R 



T\ 



II - 



+ I+- -N H (N h + 1) (5.65) 



and, observing that for the monotoniticity of the Z($) one has to choose 
linearly increasing integers as quantization numbers, hence I H = Ylj=i If = 
\Nh {Nh + 1), it follows: A = A m . This means that the generic m holes 
NLIE state characterized by minimal quantum numbers, i.e. the m solitons 
sine-Gordon state, corresponds to the primary state of the [</> m ] conformal 
family in the c = 1 BCFT. 

About the topological charge 

In formula ( J3.12I ) the topological charge of any DSG state has been intro- 
duced. We set 



^ 2tt 



.) + m + — m_ 



m± G Z 



(5.66) 



i.e. the topological charge consists of a combination of a term depending on 
boundary and an integer. It is evident that as far as the formula ( J5.591 ) is 
concerned, one could expect to give an expression of conformal dimensions 
in terms of the topological charge of the corresponding DSG state. Identi- 
fying the topological charge of the vacuum with Q = £- (0 + — <p_) and the 
topological charge of a soliton as Q\ — m — 1, the formula ( J5.64) for the 
solitonic conformal dimension becomes 



A 2ff 



Q0+2 



(5.67) 



and, in general, for any excited solitonic state over the vacuum with topo- 
logical charge Q s = m = 2S 






m 



Qo + j 



N 



(5.68) 



This means that, once the topological charges fixed, one has directly the 
corresponding conformal state in the UV limit. In particular, for any anti- 
solitonic state the topological charge comes from applying the symmetry <fi — > 
—0 and <p± — > —<p± simultaneously. Therefore, one can generalize the Eq. 
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( 15.681 ) for all antisolitonic states, which have been defined in the scattering 
theory as pure holes states under the symmetry C : h± — > — h± , S — > — S, 

as 



A = ^ 
P 2 



-Qo + % 



+ N (5.69) 



where Q_ s = — Q s = — m. Of course, due to the structure of conformal 
dimensions expression, solitonic and antisolitonic states have the same di- 
mensions. This is not surprising, in fact, the symmetries of the model makes 
irrelevant the signs of the boundary fields and of the total spin: any state is 
equivalent to the C-transformed state. This property has to be encoded also 
in the conformal limit. 

Boundary bound states as conformal states 

Among all other possible states the most interesting ones are those cor- 
responding to the boundary bound states described in the scattering the- 
ory. In order to have such boundary states one has to choose \h±\ > h c , 
i.e. the DSG boundary fields have to be in the ranges | < 0_ < g - , 

g(g±l) / X ^ ff(2p+l) 



< (f) + < i -^-| — -, according to the relations ( J5.8I ). In these ranges 
of definition, the conformal dimension of the vacuum state turns out to be 
A = ^ [(f) + — 0_] = A , where A m refers to the conformal dimension of the 
primary state in the represenatation 7i m . This conformal dimension has the 
same structure of dimension ( 15.631 ), but it is quantitatively different due to 
different values of <f>±. At this point one may distinguish several situations: 

• the case without boundary bound states can be described as before for 
all allowed states, the only difference is the value of boundary fields 
and, hence, the value of vacuum conformal dimension; 

• the case with boundary bound states which implicitely contains a vari- 
ation of boundary conditions. The description of such states is the 
following. 

The first excited state over the vacuum is described, in the large / NLIE, by 
an imaginary hole at position ia = —i^H_ (see Eqs. ( 15.171 ) and ( 15.231 )) 

with quantum number 7 (_> = and energy E = E^ s + M. cosan (notice 
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that in the UV limit the constant term E^ can be ignored because of Eq. 
( 15.421 )). It corresponds to a conformal state with Nj = 1 and conformal 
dimension 



2 



2 



(5.70) 



A soliton in the IR limit is described by a real hole $i state (iV# = 1) with 
energy Ai cosli^i over the first excited state, i.e. over the boundary bound 
state energy. This means that a real hole is no more the first excited state, 
in fact it has to describe a soliton scattering on an excited boundary, i.e. a 
free particle in a system with changed boundary states. Taking a boundary 
bound state means fixing new physical boundary states, any particle has to 
be an excitation over this reference state. The conformal state corresponding 
to a soliton in an excited boundary system has S + = S 
then 



\ and 1^ = 1, 



A 



R 



A, 



(5.71) 



In general, the Nh solitons state in a boundary bound state system corre- 
sponds to the conformal state with conformal dimensions 



1 2 



+ -(N H + Nj)R 



A 



N H +N! 



(5.72) 



As a consequence of the boundary bound states the conformal dimensions of 
the one soliton state appear to be in the Heisenberg representation Hi+at, in- 
stead of Hi. This is not surprising because a NLIE boundary bound state has 
to correspond to a physically relevant change of boundary state conditions 
and then to a change of the partition function and of the conformal theory. 
In particular, an expression for the conformal partition function taking into 
account different boundary state conditions might be 



l \ 4>+-4>- 
r2 I v^F 



■i(2S+AT/)ie 



Z DD oc -^ > VL^ 7 ^ 1 ^™^ (5-73) 

7](q) ^ 

where <p± together with Ni takes into account the boundary contributions, 
in the sense that not only the value of boundary fields but also the specific 
boundary states, i.e. the number of boundary bound states in the NLIE, 
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define the various Heisenberg representations which the conformal states are 
organized into and the partition function of the theory. It is evident that 
there are as many partition functions, therefore distinct conformal theories, 
as boundary states. Due to the NLIE, one can classify these conformal 
theories with the number Nj of boundary bound states. In formula ( 15.731 ) 
the number of boundary bound states Nj has to be fixed a priori, it depends, 
actually, on the fact that the boundary state has to be fixed before the study 
and classification of the conformal states in the representations 7i m . Let us 
observe, finally, that the integer n in the sum ( 11.731 ) has to be associated to 
the combination n = 2S + Nj of the third component of the spin and the 
number of boundary bound states. 

5.3 Conclusions 

In this thesis it has been shown how the Nonlinear Integral Equation deduced 
from the XXZ spin chain with boundary magnetic fields ( 13.201 ) describes the 
finite size effects of the sine-Gordon theory with Dirichlet boundary condi- 
tions ( J3.10) ). A summary of most important results is the following: 

1. we have presented a derivation of the NLIE ( 14.351 ) governing the fi- 
nite size effects for sine-Gordon theory with two boundaries, each with 
an independent Dirichlet boundary condition, as a continuum limit of 
the Bethe ansatz equations ( 13.321 ) of the alternating inhomogeneous 
XXZ spin chain with magnetic fields at the boundaries for the case of 
antiparallel fields. Excited states have been considered as well; 

2. by examining the infrared limit of the NLIE ( 15.11 ) it has been shown 
that it gives strong evidence that the underlying quantum field theory 
is actually the Dirichlet sine-Gordon model. The relation among NLIE 
boundary parameters H± and boundary fields <j)± of Dirichlet sine- 
Gordon has been fixed; 

3. the correct two particle sine-Gordon 5* matrices are reproduced ( 15.101 ) 
leading to a coherent description of NLIE states describing scattering 
solitonic and antisolitonic particles. In particular, the correct solitonic 
and antisolitonic Ghoshal-Zamolodchikov boundary reflection ampli- 
tudes (15.61) are obtained; 
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4. by analyzing several infrared states, a relation among sine-Gordon 
boundary bound states and NLIE states has been proposed. The cor- 
rect excited boundary reflection amplitudes ( 15.251 ) and ( 15.311 ) are re- 
produced as well and the scattering of particles with them is under 
control; 

5. by computing the ultraviolet conformal weights ( 15.551 ) from the NLIE 
it has been shown that the energy spectrum is consistent with the 
spectrum of a c = 1 boundary conformal field theory ( 15.621 ); 

6. conformal dimensions of solitonic and antisolitonic states and boundary 
bound states can be reproduced and it has been shown that they belong 
into the correct representations ( 15.611 ), even if the states which we have 
considered are not enough to describe the whole sine-Gordon spectrum, 
for example we have not taken into account breather states; 

7. we argued a general formula to relate NLIE states to boundary confor- 
mal states described by the partition function ( 15.731 ). 

The understanding of the Dirichlet sine-Gordon NLIE and of the finite size 
behaviour of the continuum theory is not quite complete. Indeed, some open 
questions have not been answered yet. Interesting questions need further 
investigation: 

1 . the whole set of NLIE solutions describing the finite volume behaviour 
of scaling functions is not completely understood. The main difficulty is 
the problem of counting, i.e. the complete knowledge of the structure 
of solutions that are highly dependent on the value of the coupling 
constant; 

2. there is an unresolved technical difficulty too. Apart the first simple 
states with only real and pure imaginary sources describing solitonic 
and antisolitonic states and boundary bound states, in order to study 
configurations of sources with non-zero imaginary parts, the under- 
standing of the NLIE boundary term over the fundamental branches of 
definition of functions ( 14.181 ) and ( J4.23I ) is required. This is crucial for 
a complete comprehension of other NLIE excited states as those which 
should describe breathers scattering off unexcited or excited bound- 
aries; 

101 



The continuum theory 



3. although the vacuum state NLIE obtained here reproduces some known 
facts (see [58], 66J, [67]), the excited states interpretation and their scaling 
behaviour should need to be numerically compared to results coming 
from other nonperturbative methods, namely the Truncated Conformal 
Space approach. The agreement of NLIE and TCS approach predic- 
tions should confirm the validity of the NLIE approach for the boundary 
sine-Gordon model. 



Finally, to conclude this thesis, let us mention that, although many progresses 
on the knowledge of general integrable boundary conditions for the sine- 
Gordon theory have been done [69, 70, [71] , from the point of view of the NLIE 
approach the understanding of the Dirichlet boundary conditions case should 
be a first step towards the investigation of the general boundary conditions 
case: first attempts have been recentely presented (see, for example, |72|, 73], 
74 , 75l , |76|), but a complete description has not been done yet. Nevertheless, 
we think that further investigation could throw some light on the general 
structure of IQFTs and finite size effects which, at this point, few people 
might consider as a pure academic or even theoretical game, even if theories 
such as sine-Gordon model are usually called toy models. 
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